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(N ; ABSTRACT 

> . 

We present u,b,r & i galaxy number counts and colours both from the North and 
| South Hubble Space Telescope Deep Fields and from the William Herschel Deep Field. 

r . The latter comprises a 7' x 7' area of sky reaching b ~28.5 at its deepest. Following 

Metcalfe et al. (1996) we show that simple Bruzual & Chariot evolutionary models 
I , which assume exponentially increasing star-formation rates with look-back time and 

■ qo=0.05 continue to give excellent fits to galaxy counts and colours in the deep imaging 
(?T) ' data. With qo=0.5, an extra population of 'disappearing dwarf galaxies is required 

, to fit the optical counts. 

t-H ■ We further find that the (r — i) : (b — r) colour-colour diagrams show distinctive 

features corresponding to two populations of early- and late-type galaxies which are 
well fitted by features in the Bruzual & Chariot models. The (r — i) : (b — r) data also 
suggest the existence of an intrinsically faint population of early-types at z ~ 0.1 with 
similar properties to the 'disappearing dwarf population required if qo=0.5. 
<i ■ The outstanding issue remaining for the early-type models is the dwarf-dominated 

Qh IMF which we invoke to reduce the numbers of z > 1 galaxies predicted at K < 19. 

JL , For the spiral models, the main issue is that even with the inclusion of internal dust 

j_i ■ absorption at the = 0.3 mag level, the model predicts too blue (u — b) colours for 

| late-type galaxies at z ~ 1. Despite these possible problems, we conclude that these 

■ simple models with monotonically increasing star-formation rates broadly fit the data 

to z ~ 3. 

We compare these results for the star formation rate history with those from the 
different approach of Madau et al. (1996). We conclude that when the effects of internal 
' dust absorption in spirals are taken into account the results from this latter approach 

are completely consistent with the T=9Gyr, exponentially rising star formation rate 
density out to z w 3 which fits the deepest, optical/IR galaxy count and colour data. 

When we compare the observed and predicted galaxy counts for UV dropouts in 
the range 2 S5> z ^ 3.5 from the data of Steidel et al. (1999), Madau et al. (1996) and 
new data from the Herschel and HDF-S fields, we find excellent agreement, indicating 
that the space density of galaxies may not have changed much between z — and 
z = 3 and identifying the Lyman break galaxies with the bright end of the evolved 
spiral luminosity function. Making the same comparison for B dropout galaxies in the 
range 3.5 Si z S5> 4.5 we find that the space density of intrinsically bright galaxies 
remains the same but the space density of faint galaxies drops by a factor of ~5, 
consistent with the idea that L* galaxies were already in place at z « 4 but that dwarf 
galaxies may have formed later at 3 Ss z S5 4. 
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1 INTRODUCTION 

* Present address: LAS, Traverse du Siphon, Les Trois Lues, F- In four previous papers, Jones et al. (1991) (hereafter Paper 

13102 Marseille, France I), Metcalfe et al. (1991) (hereafter Paper II), Metcalfe et al. 
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(1995a) (hereafter Paper III) and McCracken et al. (2000) 
(hereafter Paper IV), we used photographic and CCD data 
to study the form of the galaxy number-magnitude relation 
at both optical (B ~ 27.5) and infra-red (K ~ 20) wave- 
lengths. In this paper we extend our observations of the B- 
band field of paper III to even deeper magnitudes (B ~ 28 
mag) , and over a much larger area of 7' x 7' , by exposing for 
a total of 30 hours over five nights with the Tektronix CCD 
at the prime focus of the 4.2 metre William Herschel tele- 
scope (WHT). We have also acquired ~ 34 hours of exposure 
in the [/-band (to U ~ 27 mag), 8 hours in the i?-band (to 
R ~ 26.5 mag) and 5 hours in the /-band (to / ~ 25.5 mag) 
on the same field, again with the WHT. We refer to this field 
throughout as the William Herschel Deep Field (WHDF) . 

Although not competing in terms of either collecting 
area or field of view with these ground-based observations, 
the high resolution of the Hubble Space Telescope offers 
the ability to image sub-arcsecond objects to much greater 
depths than can be done from the ground. The Hubble Deep 
Field (HDF) project (Williams et al. 1996) has provided such 
a dataset in the public domain, and here we present our anal- 
ysis of the F300W, F450W, F606W and F8UW prelimi- 
nary release North (HDF-N) and South (HDF-S) fields, us- 
ing similar data reduction techniques as for the WHT data. 
Roughly speaking, these passbands correspond to U, B, R 
& / (see section 4.2). For unresolved objects the 3<r limit of 
B ~ 29.5 mag on the F450W image is considerably deeper 
than our WHT frame. However, for resolved galaxies this 
limit is much brighter. Indeed, for sources more extended 
than the WHT 'seeing' disk (~ 1.25" radius) the ground- 
based frame is slightly deeper, making it an ideal comple- 
ment to the HST images. A similar situation exists between 
the WHT u and the HDF F?>QQW frames, where, for a source 
at the WHT seeing limit, the WHT is actually over 0.5 mag 
deeper than the HDF. The same is not true for the F&Q&W 
and F814W images, which are 2-3 mag deeper (for typical 
galaxy sizes) than the deepest previous R and /-band galaxy 
counts (Smail et al. 1995, Hogg et al. 1997). 

The modelling and interpretation of HDF-N data has 
been discussed briefly elsewhere (Metcalfe et al. 1996). Here 
we discuss the data reduction for both ground- and space- 
based data, and present the number counts, colour distri- 
butions and models in greater detail. Finally we discuss the 
comparison between the models and the data. 

The optical - infra-red properties of the WHDF, based 
on infra-red imaging data to K ~ 20 taken on the United 
Kingdom Infrared Telescope, and data to K ~ 23 on one 
smaller, 1.8' x 1.8' area, are discussed in Paper IV. 

We now also have infra-red imaging data covering the 
whole of our WHDF to H ~ 22.5 mag taken with the Q 
Prime camera on the Calar Alto 3.5m telescope. This data 
will be discussed in a subsequent paper (McCracken et. al., 
in prep.). 



2 THE OBSERVATIONS 
2.1 New WHT observations 

Our data were taken during four observing runs at WHT 
prime focus; one of 5 nights in September 1994 during which 
all of the &-band data and about 3hrs of r-band data were 
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Figure 1. Filter /detector throughputs for the HDF (dashed 
lines) and WHDF (solid lines) observations. The WHDF curves 
include the effects of La Palma extinction. 



acquired, a further run of 4 nights in September 1995 dur- 
ing which the remaining r-band data and some of the it- 
band data were collected, and 3 nights in both October 
1996 and September 1997 when the rest of the -u-band ex- 
posures and all the i-band data were taken. Observing con- 
ditions were much better in 1994, when 4 out of 5 nights 
were photometric, than in 1995 or 1996, when the observa- 
tions were affected by dust and poor seeing. The 1997 run, 
although not completely photometric, had good seeing of 
» 1" FWHM. For the 1994 and 1995 runs we used the Tek- 
tronix CCD camera at WHT prime focus, giving a field size 
of 7.17' x 7.17' at a scale of 0.42" per pixel. The gain was set 
to 1.7e~ /ADU and the read-noise is about 7e~ . However, in 

1996 we switched to a Loral CCD, with 2048 x 2048 pixels 
at a scale of 0.26" per pixel and a gain of 1.22e~ / ADU . In 

1997 we again used a Loral CCD (although not the same 
chip as in 1996); this time the gain was set to le~ /ADU. 
These chips have a much higher (about a factor of three) 
short wavelength response than the Tek CCD, as a result 
of which most of our M-band signal comes from the 1996/7 
data. 

The Harris B, R and / filters were used - these are much 
closer to the standard photoelectric bands (throughout, any 
references to standard photoelectric R and / magnitudes 
imply those on the Kron-Cousins system) than the KPNO 
filters used in papers II and III. The U filter was a 50mm 
RGO filter, which caused some vignetting at the edges of the 
CCD. The filter response curves can be seen in Fig. To 
avoid confusion we use lower-case u, b, r & i to refer to the 
passbands defined by the WHT filter/detector combination. 

Exposures were centred on 00 h 19 m 59 s .6 +00°04'18" 
(B1950.0). This field encompasses the 24 hour INT and 10 
hour WHT images from paper III, although the centre is 
about 1.5' north of that of these previous exposures. Indi- 
vidual u and 6-band exposures were usually 2000s, with some 
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Table 1. Details of the WHT and HDF images. The WHT magnitudes are in u ccl i, fe cct j, r cc( j and i cc d 
(section 3.1); the HDF-N & S magnitudes are in the vega system (section 4.2); the STIS magnitudes are in 
the AB system. The definition of effective exposure is given in section 2.1. 
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a Inside an aperture with radius equal to the minimum radius given in Table ^ (magnitude is total magnitude 
of an unresolved object) 
b Inside 1 arcsec 2 



reduced to 1000s on occasions when the seeing was variable. 
For the r band we used an exposure time of 500s, and for the 
i band between 300 and 500s. These short exposures were 
necessary in order to avoid saturating the CCDs with the 
background sky signal. To avoid excessive atmospheric ex- 
tinction, no w-band frames were taken at an airmass greater 
than ~ 1.5, and no fe-band frames at greater than ~ 1.8. The 
r and z-band frames, where extinction is less of a problem, 
were mostly taken at higher airmass than this. The field was 
offset a few arcseconds in different directions between nights 
in order to avoid dud pixels always falling on the same place 
on the sky. The total exposure times were 59.5ks (Tek) + 
88ks (Loral) in it, 112ks in b, 27.5ks in r and 21.5ks in i. The 
effective exposure times, defined as the equivalent number 
of seconds at the zenith in photometric conditions (with the 
Loral CCD for the u) , are given in Table |j| The seeing was 
good in 1994 and 1997, always less than 2" FWHM, and 
around 0.9" on some nights, but ranged from 1" — 5" in 
1995 and 1.3" - 2.0" in 1996. The images with the worst 
seeing were not used in the data reduction (see section 3.1). 
Standard stars were taken at the beginning and end of each 
night, as were twilight flat-fields. Standards were also occa- 
sionally taken during the night, although we were also able 
to use our image frames to monitor atmospheric conditions. 

2.2 The Hubble Deep Fields 

For the HDF-N we use the F300W, FA50W, F606W and 
F814W preliminary 'drizzled' WFPC2 images released into 
the public domain on 1996 January 15 (Williams et al. 1996). 
These are flat-fielded, cosmic-ray rejected, sky-subtracted 
stacked images, resampled to a scale of 0.04"/pixel. Each of 
the three wide-field WFPC2 chips were analysed separately. 
We made no attempt to use the data from the PC chip. 
Note that some 25% of the observations are not included in 
the F450W images and an optimal weighting scheme was 
not used for stacking these data, nor for the F300W data. 
For the HDF-S we use the combined WFPC2 images re- 
leased in November 1998. We have also analysed the un- 



filtered STIS QSO field image. For unresolved objects the 
effective FWHM is ~ 0.15" (see Williams et al. 1996). 



3 WILLIAM HERSCHEL DEEP FIELD 
3.1 Data reduction 

A detailed description of our data reduction procedures can 
be found in Paper III. Briefly, a master bias frame is formed 
from a median of many individual bias exposures and nor- 
malised to the bias level in the overscan region of each ex- 
posure. This bias frame is then subtracted. The resulting 
frames are then trimmed to remove the overscan and divided 
by a master flat-field, formed from the median of all the flat 
fields in the relevant filter throughout the run. Cosmic rays 
are removed by comparing each pixel in each frame with the 
mean and standard deviation in that pixel over all other ex- 
posures (where necessary, normalised and re-scaled) taken 
that night. Any pixels more than 4.5er away from that mean 
are replaced by the mean of the surrounding 5x5 pixels 
in the image being cleaned. Finally, each image is inspected 
visually on the computer and any remaining blemishes (e.g. 
satellite trails) removed interactively. All the individual im- 
ages are then registered with one another to the nearest 
integer pixel shift by doing a pixel-pixel cross-correlation 
between images. They are then added together. Frames on 
which the seeing was equal to or greater than 2.0" FWHM 
were discarded. 

The final 6-band image is 998 x 991 pixels (7.0' x 6.9') 
and has a sky surface brightness of 21.9 mag/arcsec 2 . The 
seeing is ~ 1.25" FWHM (as in papers II and III this is 
judged by comparison with a simulated Moffat profile (Mof- 
fat 1969)). 

The final r-band image has 986 x 966 pixels (6.9' x 6.8') 
and has a sky brightness of 20.0 mag/arcsec 2 . There was 
a slight rotational shift between the 1994 and 1995 data 
which was corrected for using the shear technique described 
in paper III. The seeing on this frame is ~ 1.5" FWHM. 
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The 1995, 1996 and 1997 it-band data were reduced sep- 
arately. The 1995 data were then resampled to 0.26"/pixel 
and, along with the 1996 data, shear-rotated slightly in or- 
der to match the 1997 data. There was a problem with the 
Loral CCD used in 1997 in that it had a ~ 5% non-linear re- 
sponse at low signal levels. This posed a particular problem 
for u-band deep field data frames, where even in 2000 sec 
exposures the sky levels were low. Extensive exposure tests 
were carried out with a dome lamp to enable a correction 
factor as a function of count rate to be determined. This was 
then applied to the data. As a check on this procedure one 
standard star field was repeated with a series of exposure 
times, and detailed comparisons made between stars cover- 
ing a range of magnitudes on our field in 1997 and 1995/6. 
In fact, the uncorrected and corrected standards gave the 
same zero-point to within 0.01 mag (as expected, as most of 
the light from these stars comes at high count rates where 
the non-linear correction is negligible). However, the correc- 
tion did improve the relation between the 1995/6 and 1997 
data frames, removing a ~ 0.05 mag non-linearity between 
bright and faint stars. The three datasets were then stacked 
together (with appropriate weighting to allow for the differ- 
ent gain factors, signal strength and sky levels) and all the 
image analysis was done on this final combined frame, which 
encompasses 1593 x 1582 pixels (6.9' x 6.9') with a 'seeing' of 
« 1.35" FWHM. The sky brightness was 21.5 mag/arcsec 2 . 

The 1996 and 1997 i-band data were also reduced sep- 
arately. The non-linearity in the 1997 Loral CCD was not 
a problem for the i data, as the sky levels were well into 
the linear regime (and as we found for u data the zero-point 
from the standards was not affected). As a result no correc- 
tion was applied. However, an extra problem with all the 
i-band data was the presence of fringing. To remove this it 
was necessary to create a master fringe frame by taking the 
median of all the data-frames. This was then scaled to suit 
each individual frame (by trial and error) and subtracted. 
A small amount of rotation and pixel rescaling (~1%) were 
required before the two runs could be stacked together. The 
final i-band image has 1622 x 1747 pixels (7.0' x 7.6') and 
a sky brightness of 18.9 mag/arcsec 2 . It has the best seeing 
of all our data, with a FWHM of only 1.2". 

Both the present fo-band data and our 26hr INT prime 
focus CCD data from paper III cover the ~ 2' diameter 
area observed for lOhrs with the WHT auxiliary focus CCD 
camera (paper III). By stacking all three together in the 
common area we have been able to create an ultra-deep im- 
age. To achieve this the INT and WHT auxiliary images 
were rotated, again using the shear technique described in 
paper III, and resampled to the same pixel scale as the WHT 
prime image. The three were then registered to one pixel ac- 
curacy and added. The resulting image is 243 x 240 pixels 
(1.7' x 1.7') (as the WHT auxiliary image was circular, the 
effective exposure time is lower in the corners). The ratios 
of number of electrons detected in each separate image from 
an object are; 3.06 : 1.03 : 1.0 (WHT prime : WHT aux. 
: INT prime). The stacked image is therefore equivalent to 
1.7x the WHT prime exposure, i.e. ~ 46hrs of 4-m time, 
and is ~ 0.3 mag deeper. Fig. ^| shows a sky-subtracted sur- 
face brightness contour plot of this co-added image, with a 
geometric progression of levels starting at 1000ADU (29.4 
mag/arcsec 2 ) and rising by factors of two. 

Although not adding to the discussion of counts and 
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Figure 2. Isophotal contour plot for the 50 hour co-added 6- 
band data. Contours start at 29.4 mag/arcsec 2 and rise in 0.75 
mag steps. 



colours in this paper, this small field has the advantage of 
having been the subject of a 30 hour AT-band exposure taken 
with IRCAM3 camera (which has an almost identical field 
of view) at UKIRT. Thses data, and a detailed analysis of 
the (b—K) colour distribution, has been presented elsewhere 
(Paper IVY. 

Table [U lists the parameters of all our final images. 
3.2 Calibration 

Zero-points for our data were provided by observations of 
equatorial standard star fields from Landolt (1992). Mostly 
they were observed close to the meridian. The CCD field of 
view is large enough to see 5 ~ 10 standards on a Landolt 
field in a single exposure. 

As all our WHDF frames are on the same field we adopt 
the procedure of calibrating all the frames (in each run) in 
one band relative to ones taken on the meridian on pho- 
tometric nights. This is done by comparing magnitudes for 
a selection of bright (but not saturated) objects. The stan- 
dard stars observations then provide an absolute calibration. 
As the airmass of the standards and of the WHDF on the 
meridian are very similar, this way we are not sensitive to 
airmass extinction coefficients. In fact, by fitting a sec z law 
to the relative magnitudes of objects on our WHDF field 
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throughout a night's observations, we were able to derive 
atmospheric extinction coefficients appropriate to each run 
and hence apply relative airmass corrections to those stan- 
dards not observed at the meridian. 

Where possible, colour equations measured during the 
run are used to convert the standards to the 'natural' CCD 
bands (which agree with standard UBRI for objects with 
zero colour) and the magnitudes and colours we measure for 
galaxies are in these 'natural' systems, which we designate 
Uccd, b cc d, Tccd and i cc d (although, as we shall see below, 
these are very close to the standard photoelectric bands). 

The b cc d and r cc d magnitudes from the Tek CCD used 
in 1994 suffered a minor problem in that our observed stan- 
dard star magnitudes appeared to be correlated with expo- 
sure time (at the level of a few hundredths of a magnitude). 
After consultation with the RGO (T. Bridges, priv. comm.) 
it became apparent that a small additive shutter timing cor- 
rection was necessary for the WHT prime focus camera. The 
magnitude of this correction was dependent on zenith angle, 
being ~ +0.05s for most of our standards. Applying this 
significantly reduced the scatter between the observed and 
catalogued magnitudes, and the following colour equations 
were determined from n observations of 28 stars: 



i 1 1 1 1 i 1 1 1 

U WHT 1995 



-10 1 3 3 

(U-B) 



I 1 1 1 1 I 1 1 1 
R WHT 1995 



I 

12 3-1 

(B-R) 




Si 9 



I W11T 199? 



-10 12 3 

(R-I) 



bccd = B - 0.012(5 — R) (n = 73) 

r C cd = R [n = 65) 

The rms scatter about the B and R relations was 
±0.014 mag. The range in colour covered was —0.5 < 
(B — R) < 2.5. Fig |^ shows examples of the data on which 
these relations were based. 

In 1995 only one night was demonstrably photometric 
and this was used to calibrate our data independently from 
the 1994 observations. Again we found r cc d = R, with a scat- 
ter of ±0.009 mag from measurements of 23 standards. The 
r C cd zero-point deduced from this gave magnitudes on the 
WHDF which agreed to within 0.02 mag in the mean from 
those found in 1994. Note that no shutter correction was ap- 
plied to this data, as a series of test exposures ranging from 
Is to 20s were taken sequentially on the same field during 
photometric conditions and no dependence of the zero-point 
on exposure time was present down to the 0.01 mag level. 

The data from the u cc d and i cc d standards (whether 
taken with Tek or Loral CCD between 1995 and 1997) show 
a significantly higher scatter than the b cc d or r cc d data. The 
reasons for this are not clear; airmass corrections are not 
the cause because, as already noted, all the standards were 
taken at very similar airmass - also the effect exists between 
standards on the same CCD image and is not due to varia- 
tions from night to night (or within a single night). 

All the 1995 u-band data were tied to the one photo- 
metric night, although frequent standards were taken every 
night (the maximum variation in zero-point from night-to- 
night was 0.2 mag), and we were able to utilise most of these 
to derive an approximate colour equation of 

u T ek » U - 0.1(U - B) (n = 96) 

valid for —1.5 < (U — B) < 2.5, with a scatter of ±0.08 mag, 
excluding two standards whose magnitudes were confirmed 



Figure 3. Example standard star calibration plots (catalogue 
colour versus observed - catalogue magnitude) for the WHT u, 
b, r and i bands. The dotted lines indicate the adopted colour 
equations. 



by repeat measurements to be ~ 0.4 mag off the mean line. 
Fig ^ shows these data. 

The 1996 and 1997 u-band standards taken with the Lo- 
ral CCD showed a similar scatter with respect to the Landolt 
magnitudes of ~ ±0.10 mag. This excludes the two stan- 
dards which were discrepant in the 1995 Tek data, which 
were again found to be off by ~ 0.4 mag. Significantly, if 
we compare frames with several standards on taken both 
in 1996 and 1997, the scatter between the 1996 and 1997 
data was only ±0.03 mag. A similar exercise between fields 
in common to 1995 and 1996/7 datasets gave ±0.06 mag. 
The Loral data do not show a colour term as large as the 
0.1(1/ — B) found for the Tek data, and the scatter prevents 
any reliable estimate of a shallower relation. However, the 
standards on fields in common between the Tek and Loral 
observations indicate that there is a significant colour term 
between these two CCDs. We estimate this as 

u T ck - u Loral « -0.07(f7 - B) (n = 25) 

which reduces the scatter to ±0.04 mag. As the 1996/7 data 
dominate the stacked image in terms of signal-to-noise we 
adopt a colour equation of 

uccd ~u- om(u - b) 

as appropriate for our final dataset (for —1.5 > (U — B) > 
2.5). 

Although there is clearly some uncertainty in the u cc d 
zero-point, when the 1995, 1996 and 1997 WHDF data were 
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calibrated and reduced independently, the magnitude scales 
agreed to within 0.03 mag. 

The scatter in the i cc( j-band standards (±0.10 mag) 
appears partly due to non-photometric variations on some 
nights (especially in 1997) at the 0.05 mag level. There is also 
some evidence of inaccuracies in some of the fainter Landolt 
standards on the fields we used (and many of the brighter 
standards were unusable due to saturation), as standards on 
one field taken in both 1996 and 1997 show the same offsets 
with respect to the Landolt magnitudes (up to 0.2 mag, with 
no correlation with colour), and a scatter of < 0.01 mag be- 
tween the two years. Fig. H shows the data from 1997. As a 
result of these uncertainties, and the short colour range in 
(V — I), we make no attempt to fit a colour equation but 
just assume i cct j — I. 

Fortunately, in 1995 we observed an area encompass- 
ing the WHDF in the /-band with a Tek CCD on the INT 
2.5m telescope. These independent data show an offset of 
only ~ 0.04 mag with the WHT data, giving us confidence 
in our zero-point (although any possible problems with the 
standard magnitudes would affect these data as well). 

3.3 Image analysis 

Our procedure is very similar to that used in Papers FIV 
(note that the data in each passband are reduced indepen- 
dently); the sky background is removed and images are de- 
tected isophotally using the limits given in Table ^. These 
images are then removed from the frame, replaced by a local 
sky value, and the resulting frame smoothed heavily before 
being subtracted from the original. This produces a very 
flat background. The isophotal detection is then repeated. 
In order to reduce false detections, images whose centres lie 
only two pixels apart are recombined into one single image. 
A Kron-type pseudo-total magnitude is then calculated for 
each image, using a local value of sky. 

One complication with both our it-band and 6-band ex- 
posures was the presence of scattered light (at the 1% level 
in b and 10% in u) on one side of the field. This had a sharp 
edge which proved difficult to remove by the above back- 
ground flattening technique. It therefore proved necessary 
to remove this interactively from the frame. 

Table |^ shows the limiting isophotes and magnitudes 
for our isophotal detection routine and the minimum radius, 
radius multiplying factor (see paper III) and correction to 
total magnitude for our Kron magnitudes. As in our pre- 
vious papers the minimum radius is set to be that for an 
unresolved image of high signal-to-noise, and the correction 
to total is the light outside this minimum radius for such 
an image. Our measurement limits (table give the total 
magnitudes of unresolved objects which are a 3er detection 
inside the minimum radius. 

As in paper III we measure fixed aperture colours. For 
all the WHT data we use an aperture of 1.5" radius, and 
correct the measured colours for the difference in 'seeing' 
between the different bands (this is a small effect, and is es- 
timated from inspection of stellar profiles). Images detected 
in one band are matched with those detected in the other 
within ~ 1" of the position in the first band. Any multi- 
ple matches which occur are decided by visual inspection. If 
either image is below the 3a magnitude measurement limit 
for that band then no colour is measured. 
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Figure 4. (a) (R. — I) : (B — R) colour-colour diagrams for ob- 
jects identified as stars on our data (open circles) compared with 
the Landolt (1992) observations of standard stars (dots). Our 
magnitudes have been corrected to the standard passbands us- 
ing the colour equations in section 3.2; (b) as (a) but now for 
[U - B) : (B - R) 



We can compare (B — R) colours between our new WHT 
data and our INT results in Paper III on an object by object 
basis. We find a scatter in (B — R) rising from ~ ±0.1 mag 
for 22 < B < 24 mag to ~ ±0.5 mag for 25 < B < 26 
mag. The much noisier INT data will, of course, dominate 
the scatter. 



3.4 Star/galaxy separation 

Star-galaxy separation was done on the fe-band frame using 
the difference between the total magnitude and that inside 
a 1" aperture, as described in paper II. This enabled us to 
separate to b ~ 24 mag. Some additional very red stars were 
identified from the r and i frames in similar fashion. 

The colours of the stars can be used as an external, 
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Table 2. Parameters used in the WHT and HDF image analysis. As in Table g| WHT magnitudes are in 
the natural ccd system, whilst the HDF magnitudes are on the vega system. 

Frame Limiting isophote Limiting isophotal Minimum radius Kron multiplying Correction to total 

(mag/arcsec 2 ) magnitude (") factor (mag) 



Wxl 1 U 


on 7c: 


oq n 
zo.U 


1 Qfl 

l.oU 


1.4U 


n qo 
U.oz 


WHT b 


31.0 


29.0 


1.25 


1.40 


0.34 


WHT r 


29.0 


27.3 


1.4 


1.50 


0.30 


WHT i 


28.0 


26.8 


1.25 


1.44 


0.28 


Co-added b 


32.0 


30.0 


1.375 


1.50 


0.33 


HDF-N F300W a 












HDF-N F450VF 


28.5 


29.5 


0.35 


2.0 


0.11 


HDF-N F606W 


28.7 


29.7 


0.35 


2.0 


0.11 


HDF-N F814W 


27.7 


28.9 


0.35 


2.0 


0.11 


HDF-S F300W a 












HDF-S F450W 


28.5 


29.5 


0.35 


2.0 


0.11 


HDF-S F606W 


28.7 


29.7 


0.35 


2.0 


0.11 


HDF-S F8UW 


27.7 


28.9 


0.35 


2.0 


0.11 


STIS Unfilt. 


28.50 


31.0 


0.35 


2.0 


0.11 



a F450iy-band detections used for the F300VK-band - see text. 



additional check on the accuracy of our calibration. Fig. |j| 
shows the (U - B) : (B - R) and (R-I) :(B- R) colour- 
colour diagrams for the unsaturated stars on our frames. 
Also shown are the colours of the Landolt (1992) standard 
stars. The agreement between the stellar loci is at the 0.1 
mag level or better for all the colours. We note that we have 
no stars on the giant branch and that nearly all our stars are 
type G or later. There are, however, quite a sizeable num- 
ber of stars which lie well away from the expected locus, in 
the sense that for their (B — R) colour they are too blue 
in (U — B) and, sometimes, too red in (R — I). Visual in- 
spection reveals that most of these are isolated objects (so 
not errors due to confusion), and their appearance concurs 
with their automated identification as stellar. They tend to 
be amongst the fainter stars identified. These may be very 
compact galaxies, subdwarfs or QSO's. 
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3.5 Comparison with our previous published data 

In Fig. ^ we show the comparisons between our B and R INT 
data from paper III and our new data on this field. Colour 
corrections have been included to put all the data onto the 
standard photoelectric system. We find for the mean offset 
and rms scatter for over 500 objects 



B(whi) - B(int) = -0.10 ± 0.20 B < 25 

R(wht) - R{int) = -0.06 ± 0.19 R < 23.5 

It is clear that there is significant offset between the 
new and old data, in the sense that the new magnitudes are 
brighter. 

We have made a detailed comparison with our original 
short exposures on this field (paper II), both for aperture 
magnitudes and for Kron total magnitudes, in order to de- 
termine the source of this discrepancy. It would appear that 
the major problem comes from an error in the exposure time 
recorded in the header of one of the original B CCD frames 
on this field described in Paper II. This accounts for 0.06 
mag, and propagates to the Paper III data, as this was cali- 
brated using the results from Paper II. The remaining offset 
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Figure 5. A comparison between magnitudes from this work 
and those from our previous observations on part of this field: 
(a) WHT S-band against INT 24hr S-band, (b) WHT ij-band 
against INT 3hr ij-band. Colour equations have been used as 
indicated on the axes. 
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Table 3. Results of adding artificial stars to the real WHT data 
frames; (a) for the 30hr 6-band frame; (b) for the co-added 6-band 
frame; (c) for the 6hr r-band frame; (d) for the stacked u-band 
data; (e) for the stacked i-band data. 

True Magnitude Measured Magnitude Detection rate(%) 
(a) 



(b) 



(c) 



(d) 



(e) 



25.75 


25. oU ± U.lo 


85 


26.25 


11 in on 

2b. 11 ± 0.20 


79 


26.75 


ZO.o4 ± 0.2o 


70 


27.25 


27.23 ± 0.31 


59 


27.75 


27.83 ± 0.46 


43 


25.75 


zo.oy ± 0.21 


90 


26.25 


zo.ll ± U.zl 


78 


26.75 


OiZ CO _L n OA 

2b.u2 ± U.z4 


71 


27.25 


2 ( .lb ± U.oo 


63 


27.75 


2 ( . ( ( ± u.oy 


50 


28.25 


o o i c i n ,1 o 

28.16 ± 0.48 


46 


23.75 


23.67 ±0.13 


88 


24.25 


r\ a i r I r\ on 

z4.1o ± 0.20 


86 


24.75 


O A *7 A 1 n OT 

24.74 ± U.27 


78 


25.25 


io.oo ± U.41 


70 


25.75 


o p» to i n on 

25.78 ± 0.39 


60 


23.50 


23.42 ± 0.11 


99 


24.00 


23.94 ±0.10 


96 


24.50 


24.46 ±0.17 


94 


25.00 


24.96 ±0.16 


94 


25.50 


25.50 ±0.26 


93 


26.00 


26.05 ± 0.34 


86 


26.50 


26.53 ±0.41 


63 


22.75 


22.69 ±0.13 


91 


23.25 


23.14 ±0.18 


93 


23.75 


23.68 ±0.19 


86 


24.25 


24.19 ±0.29 


85 


24.75 


24.79 ±0.41 


75 


25.25 


25.37 ±0.51 


48 



in B, and that in R, appears to be a combination of un- 
certainties in the zero-pointing the Paper III data using the 
Paper II results, slight errors in the correction to total mag- 
nitudes in paper III and genuine calibration disagreements. 
All these effects are individually only a few hundredths of 
a magnitude, and within their respective error limits, but 
unfortunately appear to have all summed in the same direc- 
tion. 

From the identified sources of error, we believe that the 
magnitude scales for both the WHT and INT data of paper 
III should be brightened by 0.08 mag in B and 0.04 mag in 
R. Whenever subsequently referred to, these data have been 
corrected by this amount. As the counts in Paper II were an 
average of 12 fields, and the timing error was only present 
on the one frame, these results are unaffected at the 0.01 
mag level. 

3.6 Completeness corrections and simulations 

The deeper the ground-based number-counts are pushed 
the more important the corrections for confusion and in- 
completeness become. The magnitude limits in Table [j] are 
just theoretical calculated values for isolated images and are 
based on measured sky noise, and take no account of the 



efficiency of the detection procedure or the effects of confu- 
sion. To make a better estimate of the completeness of our 
counts we have added numerous artificial stars of various 
known magnitudes to our real data frames and subjected 
them to the normal data reduction procedure. Fig. ^| shows 
the distribution of measured magnitudes for these stars for 
the 30hr 6-band data. Table ^| gives the mean magnitudes, 
scatter, and detection rate for this data and for the stacked 
fa-band frame and the r-, u- and i-band frames. Note that 
an image is considered undetected if it is merged with an- 
other image and the combined brightness is a factor two or 
more greater than its true magnitude, or if it is not found 
within ±2 (Tek) or ±3 (Loral) pixels of its true position. 
As expected the detection rates in the real data drop as the 
magnitude becomes fainter. This is almost entirely due to 
images being merged with other, brighter images. Hence the 
higher detection rates in both the r-band and the i-band 
data, where the density of objects is lower, despite the sim- 
ilar signal-to-noise ratios. 

Although the above procedure gives a reasonable es- 
timate of the measurement of real images it tells us little 
about whether we are losing low surface brightness, extended 
galaxies, and nothing about the number of false detections 
we may pick up. Some of these may be genuine noise spikes, 
but most are actually caused by the effect of noise causing 
the software to deblend single images into multiple compo- 
nents. The best way to account for all these effects is to cre- 
ate full simulated CCD images, as described in paper III, for 
each evolutionary cosmological model of interest and run the 
image detection and analysis software on these. The results 
can then be compared with the real data. Ideally this would 
require a knowledge not only of how galaxies evolve in lumi- 
nosity and number, but also in morphology. In practice, such 
detailed information is not available; however, much can be 
inferred by treating galaxies as ideal bulges and disks and 
we have therefore re-created the two evolving _B-band simu- 
lations described in paper III, with the appropriate parame- 
ters for our 30hr WHT data (note that in order to simulate 
the noise correctly the simulations include a contribution 
from sources fainter than the measurement threshold). Fig. 
displays the true and measured counts from the two mod- 
els, together with that from our real data, both raw and cor- 
rected for the detection rates in Table [| Two points emerge 
from this. First, the corrections based simply on the detec- 
tion rates appear very close to those inferred from the full 
simulations. To demonstrate this further we show in Fig. ^ 
the results of a simulation in which we have adjusted the 
true count to be very close to the real count as implied by 
the detection rates. The measured count from this simu- 
lation is in reasonably good agreement with the raw data, 
apart from the faintest points where the model falls off faster 
than the data. We suspect this is a limitation of the models 
rather than an indication that the count slope suddenly rises 
again. Second, the true counts appear to lie somewhere be- 
tween our two models. In particular, that model favoured in 
Paper III, with the steep luminosity function slope at high 
redshift, is now seen to overpredict the number of galaxies. 

The corrections applied to the counts for the 30hr b- 
band frame are based on this simulation. For the other 
bands, i.e the co-added &-band frame, where with such a 
small area the corrections are particularly susceptible to the 
actual arrangement of images, and the u-, r- and i-band 
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Figure 6. The results of adding simulated stars of known magnitude to the real WHT 30 hour 6-band frame. Histograms of measured 
magnitude are shown for 162 simulations of stars of each of five known magnitudes (indicated by the dashed vertical lines). Non-detections 
are placed in the faintest bin. 



frames which are not quite as deep, we appeal to the agree- 
ment between the 6-band detection rate corrections and 
those from the full simulations and use the detection rate 
corrections determined for each band (Table ^|) to correct 
the counts. 

None of the above tests account for the problem of the 
halos of bright images, which tend to be broken into many 
faint images by isophotal detection algorithms. To some ex- 
tent the measurement of a local sky for the Kron magnitudes 
reduces this effect, but it still proved necessary to inspect all 
images around the bright stars and galaxies in our frames by 
eye and remove detections judged to be false. The effect this 
has on the counts is small (always < 10%), but it is much 
more important when attempting to measure the clustering 
of galaxies. 



4 THE HUBBLE DEEP FIELDS 
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4.1 WFPC2 Data Reduction and Image Analysis 

The HDF-N and HDF-S images were reduced in similar fash- 
ion, except where noted. 

In order to facilitate image detection, we rebinned 
the images into 0.08" pixels and trimmed off (HDF-N) or 
blanked-out (HDF-S) the regions of lower signal-to-noise 
around the edge of the frames caused by the dithering tech- 
nique used in the observations. We then followed a similar 
procedure to our WHT reductions. A 3-D polynomial of up 
to 3 rd order was fitted to the sky background and subtracted 
from the data. An isophotal image detection algorithm was 
then run over the data and a smoothed version of the data 
with the detected images removed was subtracted from the 
original. The detection algorithm was then re-run on this 
flat-background frame. We used the zero-points issued with 
the data-release as the initial basis for our photometry. 

Normally we would then use our isophotal detections 
as a basis for measuring total magnitude using our Kron- 
style algorithm. However, the HDF data, especially at the 
shorter wavelengths where for most galaxies we are looking 
far into the rest-UV, suffers from the 'problem' that most 
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Figure 7. (a) Closed and open circles represent raw and corrected 
(on the basis of the simulated star results) galaxy counts from the 
WHT b-band data compared with true (solid line) and measured 
(dot-dashed lines) counts from the two simulations discussed in 
Paper III. (b) As for (a) but now showing our best fit simulation, 
where the true galaxy count has been adjusted to give a measured 
count in agreement with the data. 



of the bright spiral and irregular galaxies are broken into 
numerous bright knots, each detected as separate images by 
the software. It was therefore necessary to re-assemble these 
galaxies, and we chose to do this by visual inspection on the 
frames of all pairs of objects whose centres were closer than 
0.35". On average, this resulted in ~ 100 re-assembled ob- 
jects per frame, composed of an average of 4 sub-images, out 
of about 2500 detected objects. Inevitably it is ambiguous 
as to whether some images should be re-assembled or not - 
many of them do not resemble 'normal' galaxies. In making 
our judgements we have borne in mind that an angular sep- 
aration of 0.1" never exceeds a true separation of ~ 400pc 
for c/o = 0.5 and is unlikely to reach even twice this for low 
go- It is difficult to conceive of images this close as separate 
entities. 

Once the 'broken' images have been reassembled we 
run our Kron magnitude software to determine total mag- 
nitudes. Tables ^ and |2| give details of the magnitude limits 
and isophotal and Kron parameters adopted for the HDF 
fields. With our choice of minimum radius (0.35") and mul- 
tiplying parameter (2.0) we should detect about 90% of 
the light from resolved or unresolved objects. However, this 
choice of minimum radius is not optimal for a star (for which 
a 0.15" radius would result in a 0.5 magnitude improve- 
ment in detection limit), but a compromise between signal- 
to-noise and the desire not to underestimate the magnitudes 
of faint, resolved galaxies. Many of these are likely to have 



their radii set to the minimum value just due to the noise in 
calculating the Kron radius. 

As the F300W frames are much less deep than the oth- 
ers (due to the poor UV response of the WFPC2 CCDs) it 
was decided to use the F450W image detections as input to 
the Kron magnitude routine for these frames, rather than 
those from the F300W image. 

For an unresolved image, the total magnitudes which 
give a 3<r detection inside our minimum radius of 0.35" are 
given in Table |l| We emphasise the point that these limits 
will be brighter for resolved galaxies. 

Although the PC data is included in the HDF-S frames 
(and went through the reduction procedure), we exclude this 
area from our final sample due to the much lower signal-noise 
on this chip. The final areas for the HDF-S and combined 
HDF-N data are listed in Table |lj These include a small loss 
caused by images being too close to the edge of the fields to 
be measured. 

Due to the vastly improved resolution and pixel scale 
the HDF images do not suffer from confusion losses due to 
crowding in the way the ground-based data does, and no 
correction has been applied. Incompleteness does occur at 
the faintest magnitudes due to the offset and scatter between 
our isophotal and total magnitudes, and due to isophotal 
effects. These are discussed in Section 4.3. 

As with the ground-based data we use fixed apertures 
to measure colours (section 3.3), but with a much smaller 
radius of 0.35". No relative 'seeing' correction is required. To 
match images we adopt a slightly different strategy to that 
for the WHT data, in that we use the positions of the images 
detected in one band as input to the Kron magnitude mea- 
suring routine on the frames in the other bands (with the 
one exception of the F300W band, where, as noted above, 
we already use the F450W detections). The reason for this 
is that, due to the high resolution of the HDF and the irreg- 
ular nature of many of the images, the positions of images 
detected in one band may not coincide precisely with those 
detected independently in another. One consequence of this 
is that we attempt a measurement of the colour for all the 
images detected in a particular band, irrespective of whether 
the corresponding object is below the detection limit in the 
other bands. 



4.2 WFPC2 Magnitude Systems 

To compare the HST results with ground-based data it is 
necessary to make some form of conversion from WFPC2 
magnitudes (throughout this paper we zeropoint the HDF 
data onto the vega system, defined so that an AO star has 
zero colour) into the standard U,B,REzI bands. This can 
only be done approximately, as none of the HDF filters are 
particularly close to their more standard counterparts, and 
in general exact colour equations have not been measured. 
Note that for our cosmological models we always use the 
correct filters and do not rely on these transforms. 

As a starting point we adopt the synthetic colour trans- 
forms of Holtzman et al. (1995), except for ^814 where we 
use their observed values. For F450W, ^6061^ and F814W 
(we defer discussion of the F300W band until the end of 
this section) the three Holtzman et al. (1995) equations we 
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use are 

B = F45CU sa + 0.23(B - V) - .003(5 - V) 2 

V = F606„ e9O + .254(1/ - I) + M2(V - I) 2 

I = F8U vega - .062(V - I) + .025(V - I) 2 

However, we really need these relations in terms of HDF 
colour, not (B — V) and (V — I), and to introduce a con- 
version for the R band. The / : F814i, e9a relation can be 
deduced in terms of (F606 — F814)„ e9a from the above 
equations, but to proceed further we use the following two 
approximate relations between the standard photoelectric 
bands which we have derived from the Landolt (1992) list 
of standard stars 

(V-R) ~ 0.58(5 -V) 
{V-I) ~ 1.95(V-ii) 

These relations are accurate to no better than 0.05 mag, but 
this is adequate for our purposes as this error gets multiplied 
by the colour coefficients, most of which are small. All five 
equations apply only for stars with (B — V) 1.4 and (V — 
I) <, 2.0. 

Combining all the above (approximating second order 
terms), we can deduce 

(B - V) « 0.94(F450 - FQQ&) vega 

and 

(V - I) » 1.44(F606 - F814)„ ega 

and hence, 

B » F450„e 9 a + 0.22(F450 - F606)„ e9a 
7? w F606„ e9a - 0.37(F606 - F814)„ e9a 
/ w F814„ e9a - 0.07(F606 - F814)„ e9a 

For the median colour of the HDF galaxies (with F8U veg a < 
28) these approximate to 

B as F450„ ega + 0.1 
i? ^ F606„ e9a -0.1 

I S3 -F814„ e9a 

Note that the F450I4 7 colour equation is almost identi- 
cal to that used in Papers I, II and III for the photographic 
bj band. Hence 

bj « F450„e 9 a- 

The F606W band is midway between the V and R pass- 
bands, and so has a large colour transform to either. Here 
we have chosen to convert to R, as there are very few pub- 
lished V-band counts. 

The F300W observations are more of a problem. This 
band is significantly shorter in wavelength than photoelec- 
tric U and the Holtzman et al. synthetic colour transforma- 
tion is large (varying by ^ 1 mag for —1 < (U — B) < 1) and 
may be multi-valued at (U — B) ~ 0.0. According to Holtz- 
man et al. there is also some question of the reliability of 
the WFPC synthetic transforms for such short wavelength 
filters. The best we can do in these circumstances is to take 
an approximate linear fit of the form 

U ~ F300 vega ~ 0.75((7 - B) 

which reproduces the Holtzman et al. results to within ~ 0.3 



mag for —1.5 < (U — B) < 1. Combined with our B-band 
transformation this implies 

U ~ F300„ e9a - 0.43(F300 - FA50) vega + 
0.09(F450 - F606)„ e9 a 

For galaxies with colours near the median of our distribu- 
tions this roughly translates to 

U ~ F300„ e9a + 0.4 

We would hope that our B,R and I transforms are accu- 
rate to ~ 0.1 mag, which is sufficient for the number-counts 
and colour distributions. Any comparison between HDF and 
ground-based U can only be described as approximate. 

4.3 WFPC2 Completeness 

The combination of ultra-high resolution but fairly unspec- 
tacular surface brightness limits can potentially lead to 
problems in measuring resolved images. To counter this we 
emphasise that our adopted magnitude limits are quite con- 
servative. For unresolved images the highest signal-to-noise 
would be achieved inside a radius of only 0.15". By choos- 
ing a minimum radius of 0.35" we are sacrificing ~ 0.5 mag 
in depth for such objects, in the hope of better measuring 
extended galaxies. Even so, the 3ct limit for a galaxy with a 
Kron radius of 1" will be ~ 1 mag brighter than that listed 
in Table pi As an example, Fig. Q shows the distribution of 
our measured Kron radii (containing ~ 90 per cent of the 
light) for the HDF-N F606W data. The solid line shows the 
3a magnitude limit as a function of Kron radius. Although 
at the faintest magnitudes the majority of measured galax- 
ies have Kron radii less than 0.6", there is a significant tail 
stretching out to ~ 1". 

We also have to consider the scatter between the isopho- 
tal magnitudes used in the detection routine (which were 
limited as in Table ^) and the Kron magnitudes finally 
adopted. As a result of this the Kron magnitudes always 
have a brighter completeness limit than the isophotal detec- 
tions. 

Taking these effects into account we consider the fol- 
lowing are reasonable completeness limits for our counts; 
F300„ e9O ~ 27; F450 vega ~ 28.5; F606„ ega ~ 28.5; 
F814„ ega ~ 27.5 (for HDF-N) and F300 vega ~ 26.5; 
F450„ e9a ~ 28.0; F606„ ega ~ 28.0; F814„ e9a ~ 27.0 (for 
HDF-S). We analysed the HDF-S data to a slightly brighter 
magnitude limit than the HDF-N due to the presence of ar- 
eas of lower signal to noise at the joins between the three 
three individual WFPC chips. 

As a check, we have run the simulation used to correct 
the WHT fo-band data (section 3.6) with parameters appro- 
priate to the HDF-N F450W data. This suggests that for 
F450W vega 29 we start to lose significant numbers of disk 
dominated (i.e. low surface brightness) objects, but that at 
brighter magnitudes the measured count is a good represen- 
tation of the true count. However, the size of images (which 
is dependent on cosmology) is a much more important factor 
in the HDF simulations than in those for the WHT, and it 
is not clear how accurate our simulations are in this respect; 
for now we are confident that corrections to our magnitude 
scale or counts will not exceed 20%. 

We have also run the detection/analysis procedure on 
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a negative image of the HDF-S F450W data. This en: 
us to look for spurious noise detections, which may ac 
the count. We find less than a 5% contribution from such 
images to the differential count at B ~ 28, and nearly all 
those which are found lie in the areas of reduced signal-to- 
noise around the edge of the field. We therefore conclude 
that such detections are not a significant problem. 



Figure 9. (a) The difference between our magnitudes and 
Williams et al. (1996) (on their AB system) for all objects de- 
tected by us on the HDF-N F450VF frames, (b) A comparison of 
the galaxy counts from the above frames - solid dots, this work, 
open circles, Williams et al.. 



4.4 Comparison with other reductions 

There are several other versions of galaxy counts based on 
the HDF frames. Ferguson (1998) has presented a compari- 
son of our HDF-N data with that of Williams et al. (1996) 
and Lanzetta et al. (1996). Clear differences exist between 
the datasets, all three of which used different photometry 
packages to analyse the HDF images. 

Here we compare in detail our HDF-N data with those 
of Williams et al. (1996), as they show the largest discrep- 
ancy with ours. Fig. ^| shows an object by object magnitude 
comparison between the two datasets for the band, 
and a comparison of the resulting counts. It is clear that 
there is a systematic magnitude offset and that our count 
slope is steeper than theirs, resulting in a difference of almost 
a factor two in the surface density of objects at the faintest 
magnitudes. Similar differences exist in the other 3 bands. 



As pointed out by Ferguson (1998) there are two reasons 
for this discrepancy. First, as noted above our magnitudes 
become systematically brighter as they become fainter. The 
exact offset is difficult to gauge as the magnitude difference 
distribution has broad wings and becomes skewed, but, as an 
example, the peak in {FGQGthis paper — F6Q6wtiUams et al) ~ 
0.05 at F606ab ~ 25, ~ 0.10 at F606ab ~ 27 and ~ 0.20 
at FG06ab ~ 28. The affect of this on the counts is actually 
quite small, generally < 10%. Second, we find objects which 
Williams et al. apparently do not detect at all. This ap- 
pears to account for the majority of the difference between 
the datasets. A visual inspection of these images leads to 
the conclusion that virtually all are merged into adjacent 
images in the Williams et al. data. We suspect Williams 
et al.'s claim that merging is not a significant problem is 
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Table 4. n-band differential galaxy counts from the WHT data. 
Corrected counts are based on measurement of artificial stars. 

Magnitude Raw N ga ; Corrected N ga ; 



\ a ccd) 


f npr frtimpi 




l UB s I 


20.75-21.25 


15 


1160 


1160 


91 9^ 91 7^ 


lo 


1 m n 


1 n 1 


21.75-22.25 


17 


1320 


1320 


22.25-22.75 


11 


3180 


3180 


22.75-23.25 


69 


5350 


5350 


23.25-23.75 


158 


12250 


12250 


23.75-24.25 


260 


20150 


20900 


24.25-24.75 


338 


26200 


27500 


24.75-25.25 


504 


39100 


41700 


25.25-25.75 


597 


46300 


50100 


25.75-26.25 


682 


52900 


67600 


26.25-26.75 


815 


63200 


100000 



optimistic, particularly as we find that even if we include 
all their split objects we still find a significant number of 
galaxies merged in their reductions but not in ours. As an 
example, in frame 2 in the range F606ab = 27 — 28 Williams 
et al. find 142 images with their fully-merged criteria (parent 
objects only) and 172 with no merging at all (only daughter 
objects). In the same range we find 197 images. Of these 49 
are missing from their fully-merged sample (at any magni- 
tude), but 31 are still missing from their unmerged dataset. 
Of these, visual inspection reveals 25 are still merged with 
close neighbours. 

As at some level the definition of what is one galaxy 
split into two and what is two close galaxies is subjective, 
it is not possible to say who is 'right' and who is 'wrong'. 
However, this ambiguity in the HDF counts should be borne 
in mind when comparing with theoretical models. 

4.5 HDF-S STIS Imaging Data 

The Hubble Deep Field South STIS, unfiltered, 43hr image 
is to date the deepest image ever taken, being about ~ 1 
mag deeper that the F606W images (although only cover- 
ing ~ 50" x 50"). The passband is effectively defined by the 
CCD and is ~ 4000A wide (FWHM) centred on ~ 6000A . 
We have analysed the data in similar fashion to the WFPC2 
HDF images. Details of the parameters used are give in Ta- 
ble |^. Areas around the central quasar and the bright star 
on one edge of the frame were removed (and replaced by 
zero signal) before analysis, and the area of the frame was 
adjusted to take this into account. We adopted the AB zero 
point of 26.39 given on the HDF-S WWW pages. As the 
passband is very different from conventional filters we make 
no attempt to convert these magnitudes into any other sys- 
tem. 



5 GALAXY COUNTS AND COLOURS 
5.1 U-band counts 

There are very few [/-band counts in the literature. Fig. 
[ic| shows the differential [/-band galaxy count for those 
we have been able to find, together with our WHT data 
from this paper (throughout, we assume galactic extinc- 
tion of E B -v = 0.02 for the WHDF, as in Paper III). All 



Table 5. Differential galaxy counts from the HDF-N and HDF-S 
F300W fields. 



Magnitude 




"NT 

N gal 


Magnitude 




"NT 

Ngal 


(F300„ ega ) 


(total) 


(deg z ) 


(F300„ ega ) 


(total) 


(deg 2 ) 


HDF-S 






HDF-N 






23.5-24.0 


19 


15000 


23.5-24.0 


31 


26500 


24.0-24.5 


31 


24400 


24.0-24.5 


47 


40200 


24.5-25.0 


47 


37000 


24.5-25.0 


55 


47000 


25.0-25.5 


65 


51200 


25.0-25.5 


91 


77800 


25.5-26.0 


84 


66100 


25.5-26.0 


108 


92300 


26.0-26.5 


114 


89700 


26.0-26.5 


137 


117100 


26.5-27.0 


169 


133000 


26.5-27.0 


164 


140200 


27.0-27.5 


183 


144100 


27.0-27.5 


231 


197000 



Table 6. b-band differential galaxy counts from the 30hr WHT 
frame. Corrections to the counts are based on the simulated CCD 
frame discussed in the text. 



Magnitude 


Raw N 


gal 


Corrected N, 


(Peed) 


(per frame) 


(deg" 2 ) 


(deg" 2 ) 


22.5-23.0 


•14 


3375 


3375 


23.0-23.5 


50 


3825 


3825 


23.5-24.0 


128 


9825 


9825 


24.0-24.5 


229 


17575 


17575 


24.5-25.0 


352 


27000 


27000 


25.0-25.5 


480 


36850 


40400 


25.5-26.0 


706 


54175 


61000 


26.0-26.5 


801 


61475 


73500 


26.5-27.0 


939 


72075 


95500 


27.0-27.5 


1064 


81650 


129500 


27.5-28.0 


1149 


88175 


205000 



the ground-based magnitudes have been converted to stan- 
dard photoelectric U. The agreement between the ground- 
based counts is good, apart from those of Songaila et al. 
(1990) which are somewhat low. These data were taken 
through a non-standard [/-filter, and the offset to standard 
U we have used (from Hogg et al. 1997) must be rather 
uncertain. Overall there seems little change in the slope of 
dlog(N)/dm ~ 0.4 from U ~ 18 to U ~ 25. Faintward of this 
there is some indication that the counts are becoming shal- 
lower. We also show as an inset our HDF F300 vega counts. 
The approximate equivalent [/-limit for the HDF data is 
27.5 mag. These appear to have quite a similar slope to the 
ground-based data, although the F300W band is much fur- 
ther into the ultra-violet than the ground-based U and the 
K-corrections for galaxies will be very different. In particu- 
lar, the Lyman-ct forest and the Lyman-limit start to dim 
the flux reaching the F300W band at lower redshifts than 
for the U band. 

Table ^ gives our WHT counts, whilst Table ^ lists our 
HDF counts as a function of FS00 vega - 

5.2 B-band counts 

Fig. [ll] shows a compilation of published B-band galaxy 
counts over the range 15 55 B 5s 29, including the corrected 
counts from our new WHT data (for the whole WHDF and 
for the co-added field) and the HDF F450W counts. All the 
counts have been adjusted to be on the standard photoelec- 
tric .B-band system. For our HDF counts the difference be- 
tween F450 vega and B is 0.1 mag - section 4.2. As in Paper 
III we have endeavoured to plot all the datasets only to their 
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Figure 10. C/-band differential ground-based count compilation, together with the HDF-N and HDF-S F300W counts (inset). The 
models discussed in the text are shown for comparison. Note that the models shown in the inset are calculated for the F300W filter, 
whilst those on the main figure are appropriate to tt cc( j. 



Table 7. 6-band differential galaxy counts from the co-added INT 
and WHT frames (equivalent to 46hrs WHT exposure). Correc- 
tions are based on the artificial star results. 



Magnitude 


Raw N 


gal 


Corrected N 


(b cc d) 


(per frame) 


(deg" 2 ) 


(deg" 2 ) 


25.0-25.5 


38 


55100 


55100 


25.5-26.0 


16 


66700 


74100 


26.0-26.5 


33 


47800 


61700 


26.5-27.0 


47 


68100 


95500 


27.0-27.5 


69 


100000 


158500 


27.5-28.0 


63 


91300 


182000 


28.0-28.5 


79 


114500 


251200 



respective 3cr limits. All the datasets have been corrected 
for galactic extinction. It can be seen immediately that our 
WHT and HDF counts are in good agreement, suggesting 
that the WHT confusion corrections are quite accurate. 

With the improved statistics and signal-to-noise over 
Paper III it is clear that the conclusion we drew previously 
is still supported and that the count slope has flattened at 
the faintest magnitudes. However, the slope is somewhat 
shallower that before, dlogN/drrib ~ 0.25, which, if this 
reflects the slope of the faint end of the luminosity function 



Table 8. Differential galaxy counts from the HDF-N and HDF-S 
F4501T fields. 



Magnitude 






Magnitude 






(F450 vega ) 


(total) 


(deg" 2 ) 


(F450 vega ) 


(total) 


(deg" 2 ) 


HDF-S 






HDF-N 






23.5-24.0 


14 


11200 


23.5-24.0 


18 


15400 


24.0-24.5 


22 


17600 


24.0-24.5 


24 


20500 


24.5-25.0 


38 


30400 


24.5-25.0 


34 


29100 


25.0-25.5 


52 


41600 


25.0-25.5 


64 


54700 


25.5-26.0 


79 


63200 


25.5-26.0 


87 


74400 


26.0-26.5 


90 


72000 


26.0-26.5 


129 


110300 


26.5-27.0 


122 


97600 


26.5-27.0 


1 19 


127400 


27.0-27.5 


158 


126400 


27.0-27.5 


191 


163300 


27.5-28.0 


235 


188000 


27.5-28.0 


279 


238500 


28.0-28.5 


299 


232000 


28.0-28.5 


389 


332500 


28.5-29.0 


232 


185600 


28.0-28.5 


485 


414600 



at high redshift as we suggested in Paper III, corresponds to 
a ~ —1.6. There is still no sign that the counts have stopped 
rising, even at B ~ 29 mag, and the total number of galaxies 
exceeds 1.5 x 10 6 /sq. deg. 

Our results for the WHDF, the WHT+INT co-added 
frame and for the HDF F450 VEga are given in Tables [], ^ 
and|| 
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Figure 11. B-band differential count compilation. The models discussed in the text are shown for comparison, calculated for our 
ground-based b cc d passband. 



Table 9. r-band differential galaxy counts from the 6hr WHT 
frame. Corrections are based on artificial stars. 



Magnitude 


Raw N 


gal 


Corrected N 


(Tccd) 


(per frame) 


(deg" 2 ) 


(deg^ 2 ) 


21.0-21.5 


23 


1825 


1825 


21.5-22.0 


11 


3275 


3275 


22.0-22.5 


59 


4700 


4700 


22.5-23.0 


101 


8050 


8050 


23.0-23.5 


154 


12275 


12275 


23.5-24.0 


224 


17875 


20400 


24.0-25.5 


373 


29750 


34700 


24.5-25.0 


493 


39300 


50100 


25.0-25.5 


696 


55500 


79400 


25.5-26.0 


833 


66425 


109600 



5.3 R-band counts 

The .R-band count compilation, shown in Fig. |l^, has 
changed significantly since Paper III. Apart from our new 
WHT data, the deepest ground-based J?-band counts pub- 
lished are those of Smail et al. (1995) and Hogg et al. (1997), 
based on observations in good seeing on the Keck telescope. 
However, HDF data provide the major change, extending 
about two magnitudes fainter than the ground-based data. 



Table 10. Differential galaxy counts from the HDF-N and HDF-S 
F606W fields. 



Magnitude 




N 9 ai 


Magnitude 






(F606„ eoa ) 


(total) 


(deg~ 2 ) 


(F606„ eoa ) 


(total) 


(deg" 2 ) 


HDF-S 






HDF-N 






23.5-24.0 


28 


21500 


24.0-24.5 


35 


29900 


24.0-24.5 


35 


26900 


24.0-24.5 


35 


29900 


24.5-25.0 


50 


38400 


24.5-25.0 


51 


43600 


25.0-25.5 


80 


61400 


25.0-25.5 


89 


76100 


25.5-26.0 


112 


86000 


25.5-26.0 


125 


106800 


26.0-26.5 


145 


111400 


26.0-26.5 


153 


130800 


26.5-27.0 


194 


149000 


26.5-27.0 


198 


169200 


27.0-27.5 


258 


198200 


27.0-27.5 


257 


219600 


27.5-28.0 


338 


259600 


27.5-28.0 


388 


331600 


28.0-28.5 


447 


361000 


28.0-28.5 


497 


424700 


28.5-29.0 


361 


277300 


28.5-29.0 


509 


435000 



All data are plotted on the photoelectric R system. To get 
to R the HDF F606 veg a measurements have been shifted by 
—0.1 mag as discussed in section 4.2. 

As with the B band, the counts continue to increase, 
with a slope dlog(N)/dm ~ 0.37 for 20 <J R <, 26, but 
with evidence of a change to a shallower slope faintward 
of this. This occurs where the HDF data take over from 
the ground-based data, and could be an artifact of the fact 
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Figure 12. ij-band differential count compilation, together with the HDF-S STIS count in the unfiltered passband (inset). The models 
discussed in the text are shown for comparison - for ground-based r cc d in the main figure and specifically for the STIS passband in the 
inset. 



that F606W is actually midway between V and R, and 
so may not have a slope appropriate to R. However, the 
mean (F606 — F814) colour for the HDF data is fairly con- 
stant faintward of 7? ~ 24, suggesting that any change in 
count slope between F606W and R is going to be mini- 
mal. The integral R counts reach ~2x 10 6 deg -2 . We note 
that the faint end slope is very similar to the B data, with 
d log N/dniR ~ 0.25. The count data are presented in Tables 



d log NJdr, 
§and[uj 



Also shown inset in Fig. [L2|are the HDF-S STIS counts. 
These probe even deeper than the HDF F606 counts and ap- 
pear to continue with a similar slope of ~ 0.25. Note that 
the apparent turn-over in the faintest bin is due to incom- 
pleteness. 



5.4 I-band counts 

This is the first time we have included an /-band count in 
this series of papers. The WHDF counts are detailed in Table 

11, whilst the HDF counts in F814 vega are listed in Table 

12. From section 4.2 we see that F814 vega ~ J. 
The counts are plotted in Fig. |l^. Note that as well as 

our HDF counts and the WHT /-counts, we also show in 



this diagram for the first time two sets of ground-based I- 
band CCD counts taken at prime focus of the Isaac Newton 
telescope. The first set, to I ~ 23, is based on 10 of the 
fields in Paper II (total area 6.3 x 10 -3 sq. deg.). These data 
consisted of 4 x 400s exposures on each field taken with an 
RCA CCD. The data reduction and analysis of these frames 
were identical to those for the B and R frames in Paper II. 
The second set are taken off a 1.5hr total exposure with a 
Tek 1024 x 1024 CCD (total area 2.4 x 10" 2 sq. deg.) centred 
on the WHDF (see also section 3.2). These reach I ~ 23.75. 
The data reduction and analysis are similar to that described 
here for the WHT data. Once again, all magnitudes have 
been adjusted to be on the standard ground-based system. 

The ground-based counts have a slope of dlog(N) / 'dm ~ 
0.33 for 21 ^ / ^ 25, which appears to flatten slightly to 
dlog(N)/dm ~ 0.27 for the fainter HDF data. 



5.5 Colours 

As noted in sections 3.3 and 4.1 we use fixed aperture 
colours. However, the difference in this aperture between 
the WHT and HDF data (1.5" radius as opposed to 0.35") 
certainly means that we are measuring colours (on average) 
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Figure 13. I-band differential count compilation. Again, the models discussed in the text are shown, calculated for our ground-based 
i ccd band. 



Table 11. i-band differential galaxy counts from the WHT data. 
Corrections are based on artificial stars. 



Magnitude 


Raw N 


gal 


Corrected 


(iced) 


(per frame) 


(deg" 2 ) 


(deg" 2 ) 


20.0-20.5 


17 


1175 


1175 


20.5-21.0 


18 


3320 


3320 


21.0-21.5 


5(5 


3870 


3870 


21.5-22.0 


79 


5460 


5460 


22.0-22.5 


125 


8640 


8640 


22.5-23.0 


205 


14170 


15500 


23.0-23.5 


283 


19560 


21900 


23.5-24.0 


413 


28550 


34700 


24.0-25.5 


490 


33870 


42700 


24.5-25.0 


623 


43070 


61700 


25.0-25.5 


585 


40440 


85100 



over a different metric aperture from the two data-sets. This 
should be borne in mind when interpreting the data. 

Fig ^ displays the b cc a versus (b — r) C cd colour magni- 
tude histograms for our whole 7' x 7' field, split into four 
magnitude ranges. Fig |l5| shows the equivalent F^50 vega - 
limited (F 450 - F606) vega data for the HDF-N and HDF-S. 
The percentage colour completeness is indicated on each of 
the histograms. This is particularly severe in the faintest 



Table 12. Differential galaxy counts from the HDF-N and HDF-S 
F814VK fields. 



Magnitude 






Magnitude 






(F814„e ff o) 


(total) 


(deg" 2 ) 


(F814„ e9 a) 


(total) 


(deg" 2 ) 


HDF-S 






HDF-N 






22.5-23.0 


31 


23900 


22.5-23.0 


22 


18800 


23.0-23.5 


27 


20800 


23.0-23.5 


30 


25600 


23.5-24.0 


39 


30100 


23.5-24.0 


15 


38500 


24.0-24.5 


50 


38600 


24.0-24.5 


62 


53000 


24.5-25.0 


83 


64000 


24.5-25.0 


92 


78600 


25.0-25.5 


117 


90200 


25.0-25.5 


136 


116200 


25.5-26.0 


161 


124100 


25.5-26.0 


169 


144400 


26.0-26.5 


212 


163500 


26.0-26.5 


214 


182900 


26.5-27.0 


255 


196600 


26.5-27.0 


282 


241000 


27.0-27.5 


315 


242900 


27.0-27.5 


437 


373400 


27.5-28.0 


335 


258300 


27.5-28.0 


547 


467400 


28.0-28.5 


194 


149600 


28.0-28.5 


573 


489600 



ground-based bin due to the comparatively bright limit of 
the r cc d data. Also shown are the predicted histograms 
for two of our evolutionary models (these are discussed in 
sections 6.3 and 6.4). Figs and [j^ show similar plots, 
but now for b cc[ j-limited (it — b) cc d and i^SO^eoa-limited 
(F300F450)„e 9 a colours. Figs [jj and |l| show the plots for 
(r — i) C cd and (F606 — F8U) veg a colours. 
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The colour-colour diagrams for the WHDF, (r — i) C cd '■ 
(b — r) cc d and (u — b) C cd '■ (6 — r) cc d, are shown in Fig. 
BCl(a),(b). Fig. E3J shows the equivalent plots for the spaced- 
based data, split into bright and faint magnitude ranges. We 
discuss these in more detail in the next section. Note, how- 
ever, that due to the different K-corrections between the 
HST and WHT passbands, especially in the ultra-violet, 
where the filters differ by ~ 600A , the tracks of galaxies 
in these plots are not expected to be the same. 

6 GALAXY EVOLUTION MODELS 
6.1 Overview 

In our previous work (Shanks 1990, Paper II, Paper III) 
we have noted that if the B-band models are normalised at 
B ~ 1 8 rather than B ~ 15 then non-evolving models give a 
reasonable representation of the B counts and redshift dis- 
tributions in the range 18 ^ B ^ 22.5. Some support for this 
high normalisation comes from work on HST galaxy counts 
subdivided by morphology, where non-evolving models with 
high normalisation give an excellent fit to both spiral and 
early-type counts with 17 ^ I ^ 22 (Glazebrook et al. 1995, 
Driver et al. 1995). In Paper IV we have also shown that 
these same non-evolving models fit the 7f-band counts as 
faint as K ~ 23. The B counts now have the work of Bertin 
& Dennefeld (1997) added at the bright end, comprising pho- 
tographic measurements over 145 deg 2 . A gain these show a 
steeper slope than the models at B < 17. The same au- 
thors' R counts show a similar effect at R < 16. In the / 
band the recent counts of Postman et al. (1998) also show 
a steep slope at I < 17. Although these counts cover only 
16 deg 2 , they are made using a CCD. Finally, the 2MASS 
counts (Cutri & Skrutskie 1998) which cover 158 deg 2 to 
H = 15 and K = 14 are also found to be steeper than 
the no evolution models (see McCracken et. al. in prepara- 
tion). This effect therefore seems to be present over a wide 
range of passbands. If number count steepness were caused 
by evolutionary changes in star formation rate (SFR) at low 
redshifts, z ^ 0.1, then it would be surprising if the near-IR 
counts were affected as much as the B counts. We therefore 
continue to believe that the most likely interpretation is that 
this steepness at bright magnitudes is caused by large scale 
inhomogeneities in the galaxy distribution on ~ 150h _1 Mpc 
scales (Shanks et al. 1990, Paper III). Further evidence for 
this has come from preliminary luminosity function data in 
the 2dF Galaxy Redshift Survey which shows a LF with the 
local M* and slope a (Efstathiou et al. 1988, Loveday et al. 
1992, Ratcliffe et al. 1998), but a 0* which is 50% higher 
than the local value and in good agreement with the value 
that is used here. We therefore believe that these develop- 
ments justify our normalising the count models at B ~ 18, 
corresponding to an overall value of <f>* ~ 2.4 x 10 -3 . The 
exact normalisation ((f)*) we use is given in Table |l3| for the 
evolutionary models. In the case of no evolution we increase 
these <f)*'s by 10% to give the same normalisation at B — 18 
as in the evolutionary models. 

Metcalfe et al. (1996) presented two evolutionary mod- 
els with this high normalisation which gave reasonable fits 
to the counts at all the wavelengths for which they had data, 
and to the HDF colours. These were a simple, qo = 0.05 lu- 
minosity evolution model with exponentially decaying star 



formation rates, and a qo = 0.5 model, similar but for the 
addition of a population of dwarf (hereafter 'dE') galaxies 
with a constant star-formation rate until z ~ 1 and zero 
thereafter. 

As in Metcalfe et al. (1996), for simplicity we only con- 
sider two basic forms of evolution; a r=2.5Gyr exponentially 
decaying star formation model for E/S0 and Sab galaxies, 
and a r=9 Gyr model for Sbc/Scd/Sdm, both from Bruzual 
& Chariot (1993). This assumes that the evolution of the dif- 
ferent morphological types is simply governed by their SFR 
history. For qo=0.05 and qo=0.5, the present day galaxy 
ages are assumed to be 16Gyr and 12.7Gyr, implying for- 
mation redshifts of zj — 6.3 and Zf = 9.9. We also consider 
a model with a cosmological constant, zero spatial curvature 
and f2o=0.3, with a present day age of 18 Gyr and a forma- 
tion redshift of Zf = 7.9. Details of the present day lumi- 
nosity functions and colours used in our models are given in 
Table |l^. All models take account of Lyman-a forest /break 
absorption as described by Madau (1995). Similar models 
have been used by Pozzetti et al. (1998). 

One feature of these pure luminosity evolution (PLE) 
models is that they predict that z > 1 galaxies should be 
detectable at reasonably bright magnitudes and this predic- 
tion has been broadly confirmed by B ~ 24 redshift surveys 
(Cowie et al. 1995, 1996) undertaken on the Keck 10-m tele- 
scope. However, to provide a detailed fit to the B-band n(z) 
distribution it is necessary to moderate the number of high 2 
spiral galaxies predicted by the r~9 Gyr model by including 
the effect of internal dust extinction, with A\ oc 1/A (Met- 
calfe et al, 1996, cf.Campos & Shanks 1997 and refs. therein). 
We have, at least initially, assumed that the extinction law 
and the star- formation history /metallicity evolution can be 
treated independently. Edmunds & Phillipps (1997) and Ed- 
munds & Eales (1998) have investigated the effects of drop- 
ping this assumption and indeed in Section 7 of this paper 
we discuss possible evidence for evolution of the extinction 
law with redshift. A similar problem exists in the if-band 
n(z) with the r=2.5Gyr model for early-type galaxies. Here, 
essentially passive evolution overpredicts the high z tail at 
K ~ 19. We have overcome this by adopting a dwarf domi- 
nated IMF, with a slope x = 3, for early-type galaxies (see 
Paper IV). In the optical bands this produces results very 
similar to the more widely used Salpeter or Scalo IMF's. 
There is some evidence that passively evolving models over- 
predict the HST morphological counts of early-type galaxies 
(see e.g. Driver et. al. 1998 and references therein), but our 
model does a reasonable job of matching these counts, ex- 
cept at the faintest magnitudes, where morphological identi- 
fication becomes difficult. One problem for the x=3 model is 
that it predicts an M/Lb ~120M/Lq for early-type galaxies 
as opposed to M/L_b=20 for x=1.35 and compared to the 
observed value of M/L_e = lOhM/L0. This disagreement in 
the x=3 case is caused by the excess of low-mass stars. We 
have therefore also considered a model still with r=2.5Gyr 
but with an x=3 IMF which cuts off below 0.5M Q . The 
isochrone synthesis code for this case was specially supplied 
by G. Bruzual. This model then predicts M/L_b=5M/L0, in 
much better agreement with observation. The effect on the 
galaxy colour predictions is reasonably small; for example, 
at 16Gyr the predicted rest (V — K) colour is 3.27 for the 
x=3+0.5 M© cutoff case compared to (V — K) = 3.67 in the 
x=3 uncut case and to (V — K) = 3.25 in the Salpeter case. 
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Values observed for giant ellipticals are in the range 3.1-3.5. 
Vazdekis et al. (1997) have further shown that dwarf dom- 
inated IMF's with a low mass-cut, solar metallicities and a 
16Gyr age do produce consistent population synthesis model 
fits to nearby early-type galaxy spectra. Another problem 
for a dwarf-dominated IMF is that it contains too few gi- 
ant stars to generate the observed metallicity of present-day 
early- type galaxies. However, this is also a problem for a 
Salpeter IMF, and a short initial burst with a much flatter 
IMF slope usually has to be invoked to solve this problem, 
(e.g Vazdekis et al. 1997). Indeed, Vazdekis et al. find that 
their chemical evolution, 'closed-box' model with a flat IMF 
during an initial short burst and then a steeper IMF(x=2- 
3) plus low mass cut-off, produces an excellent fit to the 
spectra, and hence the observed metallicities, of local giant 
ellipticals. 



6.2 Galaxy Count Models 

Galaxy counts generally constrain particular combinations 
of evolutionary models and cosmological parameters. Figs 
[l^ - show how our various models fit to the counts in the 
four bands. The HDF-S STIS count models (see Shanks et 
al. in prep.) are shown inset into Fig. |l2|. As well as the two 
evolving models described above, we also show their non- 
evolving counterparts, with the exception of the A model 
which is very similar to the qo=0.05 no-evolution model. 
Also shown is the go = 0.5 evolutionary model without the 
extra dE component and a new qo=0.5 evolutionary model 
which has a steeper luminosity function slope for Scd/Sdm 
galaxies (see below). As can be seen, even at the depths 
reached by STIS, our preferred evolutionary models fit the 
counts reasonably well over a wide range of wavelengths. 
This agreement continues into the infra-red (Paper IV). 

Note that both the evolution and no-evolution models in 
the range 18;$ U $24 have different slopes when expressed 
in the ground-based U and HST F300W passbands, with 
the slope in the F300W band being steeper. This is due 
to the HST band peaking at 3000A and the ground-based 
band peaking at 3800A . 

As the counts move bluewards from I to B, the pat- 
tern is that the evolutionary excess of galaxies sets in at 
increasingly lower space densities of galaxies, indicating the 
increased sensitivity of the bluer bands to the increased SFR. 
However, at U the evolutionary excess appears smaller than 
at B; it also sets in no earlier in terms of galaxy number 
density. The reason is that the U band is unique in that it 
lies shortward of the 4000 A break even at z — 0. Thus the 
K-corrections for late types are less than in B by about a 
magnitude at z > 0.5 and indeed are negative for z > 1. 
This means that the no-evolution U models are intrinsically 
steeper than in B at U ~ 24 and lie closer to the U data, 
which in fact has a similar slope to the B data. Thus less 
U evolution might, in principle, be thought to be needed. 
The Bruzual & Chariot models, however, predict a similar 
amount of evolution in U as in B, about 1 magnitude at 
z > 1. This leads to the U evolution models being slightly 
higher than the data at U ~ 24 whereas the B evolution 
models are slightly lower than the data. This results in a 
problem for the models in fitting the detailed (U — B) colours 
at B <; 24, as we discuss further below with regard to the 
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Figure 14. b cc[ j-limitod (6 — r) cc( j colour-magnitude histograms 
for the WHDF data, split into four magnitude bins. The percent- 
ages indicate the colour completeness for that magnitude range. 
The predictions of the go = 0.05 (solid line) and go = 0.5 dwarf 
(dashed line) models are also shown. 



n(u — b) cc d histograms in Fig. |16|and the (u— b) cc d '■ (b—r) cc d 
distribution in Fig. po|(b). 

The qo=0.5 evolutionary model (without the added dE 
component) fits the counts in the range 18 $ B $ 25 
whereas the qo=0.05 and the similar k — 0, Qa=0-7, f2o=0.3 
models extend the fit to 18 B $ 27. It is well known that 
the 4x smaller volume available to z = 4 in the qo=0.5 case 
makes it impossible to fit the faintest optical data with a 
PLE model (Yoshii & Takahara 1988, Guiderdoni & Rocca 
Volmerange 1990, Koo et al. 1993). Metcalfe et al. (1996) 
therefore invoked the extra dE population at high redshift 
to improve the fit of this model. The parameters for this pop- 
ulation are also given in Table Essentially, these galax- 
ies have a constant star-formation rate between Zf = 9.9 
and 2 = 1 when the SFR cuts off. This causes the galaxies 
to dim dramatically in brightness, till they end up at the 
present day with an M* some 5 magnitudes fainter than the 
normal M*(z = 0) and with the colour of an E/SO. This 
is in the spirit of a disappearing dwarf model (e.g. Babul 
& Rees 1992, Babul & Ferguson 1996, see also Phillipps & 
Driver 1995 for a discussion) built in the Bruzual & Char- 
lot 'single-burst' framework. Clearly it is not unique and if 
qo=0.5 then other models such as dynamical merging may 
also produce the increased numbers of galaxies required at 
faint magnitudes. 

A further qo=0.5 model that we consider explicitly here 
for the first time is a model where we use a steeper faint end 
slope for the Scd and Sdm luminosity functions, a — —1.75 
rather than a = —1.5. Although this slope is steeper than 
seen in local luminosity functions, if the local under-density 
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Table 13. Details of the luminosity evolution models (Hq = 50kms~ 1 Mpc~ The dE type is only employed 
in the go = 0.5 case. <j>*s are given for our evolving models; 10% higher </>*s are used for the no-evolution 
models. 
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Figure 15. As Figure but now for the HDF-N and HDF-S 
F450„e 9 a-limited (F450 - F606)„ ega data. 



is real then it could be that the local luminosity function is 
deficient in faint galaxies due to biasing. Then even cluster- 
free luminosity function estimates may give misleading re- 
sults, assuming as they do that the luminosity function form 
is independent of galaxy environment. Alternatively, the late 
type galaxy luminosity function may be assumed to evolve 
to steepen from a = —1.5 to a = —1.75 at z > 1. Possi- 
ble evidence for such steepening at lower redshift has come 
from Lilly et al. (1995). Such a model gives a good fit to 
the optical UBRI counts in Figs. |l(j - [ll|. We also find that 
it is in reasonable agreement with the galaxy n(z) data at 
B < 22.5 and 22.5 < B < 24. McCracken et. al. (in prep.) 
find that the H counts are only slightly over-predicted by 
such a model at H — 28. As they point out, the flat slope of 
the H counts seems to disallow any such steepening of the 
luminosity function in the qo=0.05 cosmologies. 



Figure 16. b cc[ j-limitcd (u — b) cct i colour-magnitude histograms 
for the WHDF data. The percentages indicate the colour com- 
pleteness for each magnitude range. The predictions of the qo = 
0.05 (solid line) and go = 0.5 dwarf (dashed line) models are also 
shown. 



6.3 Galaxy Colour Distributions 

The matching of both the qo=0.05 and the qo=0.5 (dE and 
steep LF) evolutionary models to the count data over a wide 
wavelength range does not imply that the colour distribu- 
tions predicted by the models will agree with the data, as 
the colours are more sensitive to model details. In Figs. |l6| 
- ^| we show in the form of histograms the comparison be- 
tween our models and the data in both the HDF and WHDF. 
Generally we regard the agreement in all bands as good, 
at least in terms of where the broad body of the galaxy 
colours lie. The exception is in the UV, where both the HDF 
(F300 - F450)„ ega and the WHDF (u - b) ccd colour data 
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Figure 17. As Figure |uj but now for the HDF-N and HDF-S 
F450„ e9a -limited (F300 - F450) vega data. 



Figure 19. As Figure [l| but now for the HDF-N and HDF-S 
F606„ e9a -limited (F606 - F8U) vega data. 




Figure 18. r cc[ j-limited (r — i) cct j colour-magnitude histograms 
for the WHDF data. The percentages indicate the colour com- 
pleteness for each magnitude range. The predictions of the <jrj = 
0.05 (solid line) and go = 0.5 dwarf (dashed line) models are also 
shown. 



peak some 0.5 magnitudes redder than do the models. Recall 
that the spiral types in our models are already assumed to 
have dust absorption, Ab = 0.3 mag at z — 0. If we added 
more dust in the context of these models then we would have 
difficulty in matching the extended n(z) in the Keck redshift 
surveys at B < 24 (Metcalfe et al. 1996). We shall also ar- 
gue later that simply invoking dust extinction laws with a 
2200A dust feature as seen in the Galaxy is not supported 
by the (u — b) cc d '■ (b — r) cc d colour-colour data. The only 
other parameter which can potentially affect the UV colour 
of a star-forming galaxy is the IMF. We shall look further at 
what changes might improve the fit of the (u — b) cc d colours 
in Sect 6.4. 

Metcalfe et al. (1996), Madau et al. (1996) and Shanks 
et al. (1998) all assumed on the basis of the Bruzual & Char- 
lot models that in the HDF (F300 - F450)„ ega > corre- 
sponded to z ii 2 galaxies. However, the above problem, 
together with the sensitivity of U — B colours to both dust 
and Lyman absorption encourages us to use the extensive 
spectroscopic data in the HDF to check the reliability of our 
models. We shall see in Fig. g(b) that (F300- F450) vega =0 
clearly discriminates spectroscopic galaxy redshifts, above 
and below z = 2. Fig. |B| also shows that the HDF UV 
dropout galaxies accurately fall in their predicted place for 
z £ 2 galaxies in the (F450- F606)„ ega : (F606-F814)„ e9 a 
diagrams, which strongly suggests that z <; 2 galaxies can 
be discriminated from z < 2 galaxies on the basis of their 
ground-based B — R : R — I colours (see also fig. 4a of Met- 
calfe et al. 1996). Thus it seems that, whatever the model 
uncertainties at z — 1, by z = 2 the effects of Lyman absorp- 
tion dominate, leaving the model (F300 — F450) vega colours 
roughly correct at z — 2. Metcalfe et al. then went on to de- 
termine the fraction of HDF UV dropouts at 27 < B < 28 
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and found ~47±7%, a fraction which compares well with 
that predicted by the models. The results of calculating the 
z Jj 2 UV luminosity density with these criteria are discussed 
in Section 6.5 below (and by Shanks et al. 1998). 

In the WHDF a simple (it — b) cc d cut is not adequate to 
discriminate high z galaxies, partly due to the difference in 
passbands, but mainly due to the brighter magnitude limit 
and hence lower average redshift of the sample. As a result 
it is necessary to apply a simultaneous cut in (b — r) cc d, in 
order to eliminate low z early type galaxies (especially for 
r- limited samples). In order to compare with previous work 
in the literature we take as a starting point UV dropout cri- 
teria which match those of Steidel et al. (1999) as closely 
as possible; these transform to a cut with approximately 
(u — b) cc d >1 (and —0.1 < (6 — r) cc d < 1.6) in our mag- 
nitude system, and a limit of r cc d <25.25. We also assign 
(u — b)ccd colours (which are really blue limits) to those 
galaxies which were not detected in u by assuming their 
u C cd magnitude to lie at the 3<r limit (u cc d = 26.8). Given 
the relative depths of our tt-band data and Steidel et al.'s, 
this is roughly equivalent to the la magnitude they adopted 
for their non-detections. Many of these galaxies then fail 
the selection criteria due to being too blue in (u — b) cc d- 
The main purpose of this procedure is to remove faint, low- 
z early-type galaxies which might otherwise be mistaken for 
high-z dropouts. However it should be noted that it also 
loses potential high-z candidates in the process. We shall 
see how some of these might be recovered in section 6.4. 

With these criteria (and an eyeball check to exclude 
photometry problems such as mergers) we detect 32 UV 
dropouts in 49 arcmin 2 , i.e. 0.65±0.12 arcmin -2 , to r cc d < 
25.25 which is comparable to the 1.18±0.04 arcmin -2 sky 
density for z ~ 3 galaxies found by Steidel et al. (1999) at 
Rab < 25.5 (or r cc d <25.25) and the upper limit of 3.3±0.3 
arcmin -2 for UV dropouts set by Guhathakurta et al. (1990) 
at a similar _R-band limit. Applying the same blending cor- 
rections assumed by Steidel et al. (1999) for the brighter 
galaxies would increase these three observed sky densities 
by a factor of 1.3. 

6.4 WHDF Colour-Colour Distributions 

We first discuss the WHDF (r — i) cc d : (6— r) cc< j colour-colour 
diagram. This is shown in Fig. M(a). The data shows two 
unmistakable features. The first is the plume of red galaxies 
which moves redwards from (b — r) cc d ~ 1.6 to (6 — r) cc d ~ 3 
at (r — i) C cd ~ 0.7. The second is the 'hook' feature which 
lies below this red plume and contains the majority of the 
galaxies. Interestingly, the models seem to show the same 
generic features, mostly due to the redshifting of the H/K 
break through the passbands. The red plume is identified 
with early type galaxies, with their decreasing UV spectrum 
shortward of the Ca II H/K break moving through the b 
band and making the (b — r) cc< j colour redder until z ~ 0.5. 
At 0.5 ;S z 0.9 the H/K feature moves through the r band 
making the (b — r) cc d colour more blue and the (r — i) cc d 
colour more red. This shift of the (6 — r) cc d colour to the 
blue is helped by the r=2.5 Gyr evolution whose effect on 
the UV spectrum is increasingly felt in the b band at z <; 0.5. 
The hook feature is identified with spiral galaxies, with again 
the H/K break moving through the b band until z = 0.25 
when the increasing spiral UV spectrum is redshifted into 



the b band. Then, helped by the H/K break moving through 
the r band at z i; 0.5, (b — r) C cd continues to become bluer 
until the Lyman a forest enters the b band at z ~ 2.3. The 
H/K break entering the r band at z ~ 0.5 causes (r — i) cc d to 
become redder before the same feature entering the i band 
at z > 0.8 causes (r — i) cc d to become bluer again. 

The fact that the models show these features suggest 
that they are broadly correct. Indeed the evolutionary mod- 
els show increased agreement with the data than simply us- 
ing the K-corrections, especially in the case of the early- 
types beyond z ~ 0.8. This is shown in Fig. |5l] where we 
compare the colour tracks from our evolutionary models to 
the K-corrections calculated from our Bruzual models. In 
the case of the spirals, the K-corrections provide a better 
first-order description of the colour-colour tracks than for 
the early-types but evolution at low redshift and dust red- 
dening at high redshift also help; evolution makes the spi- 
rals slightly bluer than the K-corrections in (b — r) cc d and 
(r — i)ccd at z ;S 1 but then the K-corrections and evolu- 
tionary model colours become very similar in the 1 £ z Si 3 
range. The effect of spiral dust actually makes the evolu- 
tionary model redder in (6 — r) cc d at z <; 1.5 and here it is 
the effect of evolution and dust combined which improves 
the fit of the Sbc model to the data. 

However, in detail there remain differences between the 
models and the (r — i) cc d '■ (b — r) cc d data. The most serious 
problem in Fig. EC](a) is that the early-type galaxies appear 
to be too blue in (r — i) cc d at z <; 0.5, indicating more flux 
shortwards of the H/K break than predicted in the models. 
This problem is a feature of all models tested including those 
with a Salpeter IMF. Increasing r helps by extending the 
period of active star-formation to lower redshifts but then 
the (6 — r) cc d colours tend to turn to the blue at z = 0.5 
too quickly in comparison with the data. It is possible that 
models with a less sharp H/K break at high redshift may 
help because of decreased metallicity but so far we have not 
tested this hypothesis further. 

Fig. ^ compares the predicted number density of galax- 
ies at each point in the WHDF (r — i) cc d '■ {b — r) cc d 
colour-colour diagram with the data. It can be seen that 
the qo=0.05 and qo=0.5 dE evolutionary models both gen- 
erally give a good fit. In both cases the spirals are well fitted. 
As in Fig. |20|(a) , the evolutionary models predict early-type 
colours which are too red in (r — i) C cd for z > 0.5. It should 
be noted that in the case of the early-types the qo=0.5+dE 
model appears to give more low redshift early types and 
thus better agreement with the data than the low qo model. 
These galaxies are estimated to lie at z ~ 0.15 from their 
colours and thus will have Ms ~ — 16. These low luminosity 
galaxies are therefore excellent candidates for a disappearing 
dwarf population and spectroscopy is required to confirm 
their redshifts and investigate the nature of these galaxies. 

The 32 UV dropout, r ccd < 25.25, candidate 2.5 <, z ^ 
3.5 galaxies based on the Steidel et al. (1999) selection crite- 
ria discussed in the previous section have been circled in the 
(r — i)ccd '■ {b — r) cc d plane in Fig. However, it is clear that 
we can use (r — i) cc d to refine further our candidate selec- 
tion, as some of the dropout galaxies have (r — i) cc d colours 
that are too red for them to be at high redshift. At first 
glance this would seem to reduce our numbers even further. 
However, if we re-examine those galaxies with no measur- 
able u magnitudes which were assigned u cc d magnitudes at 
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Figure 20. (a) fe cc[ j-limited (r — i) cct j : (b — r) cc[ j and (b) (u — b) ccd : (b — r) oc( j colour-colour diagrams for the WHDF data, showing the 
evolutionary tracks predicted from our Bruzual & Chariot models for the five galaxy types for < z < 3.5. Boxes indicate 0.5z intervals. 
The tracks are colour coded according to redshift; < z < 0.5, blue; 0.5 < z < 1, cyan; 1 < z < 2, green; 2 < z < 3, red; z > 3, magenta. 
Arrows show limits to colours for galaxies not detected in i, r or ti, assuming the missing object has a magnitude at the 3<r limit for each 
band. The enclosed, dotted area on the (u — b) cc( i : (b — r) ccd plot shows the u-dropout selection of Steidel et al. (1999) transformed into 
our passbands. The locations of our initial sample of 32 high-z, UV-dropout galaxies based on the criteria of Steidel et al. (1999) are 
shown on both plots as circles and inverted triangles, the latter symbol indicating candidates which pass the Steidel et al. criteria but 
are too red in (r — i) cc[ j. The extra candidates which increase our final sample to 43 dropouts are shown as upright triangles (note that 
all the candidates are selected from an r-limited sample and so may not have corresponding data points or arrows on these plots). 



the 3<r limit and subsequently failed Steidel et al.'s selection 
criteria (see section 6.3), we find that many of these have 
(r — i) C cd colours which indicate they are 2.5 z 3.5 
galaxies ((r — i) C cd < 0.2). The net result is that our number 
of dropout candidates increases from 32 to 43, or 0.88 ±0.13 
armcmin~ 2 . It is this sample that we take forward to our dis- 
cussions in Sections 6.7 and 6.8. The qo=0.05 and 0.5 evolu- 
tionary models predict a sky density of 2.6-2.9 arcmin -2 for 
2.5 Ss z 3.5 galaxies at r cc d <25.25. Thus the sky densities 
predicted by both models are in reasonable agreement with 
the sky density of z ~ 3 galaxies found in the WHDF and by 
Steidel et al. (see section 6.3), especially if we apply Steidel 
et al.'s bright-end 30% blending correction (which should be 
relevant to both datasets). This suggests that Lyman break 
galaxies may be interpreted as the high luminosity tail of 
normal spiral galaxies as their SFR increases exponentially 
back towards their formation epoch. This conclusion applies 
both in the qo=0.05 case where 99% of z <; 3 galaxies are 
Sbc/Scd/Sdm and even in the qo=0.5 'disappearing dwarf 
dE case where 86% of the z <; 3 galaxies are Sbc/Scd/Sdm 
and the remainder mostly E/SO/Sab with less than a 1% 
contribution from dE. 

It can also be seen from Fig. |2^(a) that many other 
galaxies appear in the same (r — i) C cd '■ (b — r) cc d location 
as our UV-dropout candidates but with no UV dropout. 
Clearly there is the possibility that these are also z ~ 3 
candidates. Madau (1995) has used the incidence of Lyman 



limit systems in QSO spectra to argue that almost all (~ 
80%) z <; 3 galaxies should show UV dropout. This argu- 
ment works well at z ^ 3.5, because there is a relatively 
large redshift range between 3 ^ z ^ 3.5 for a Lyman limit 
system to enter the U band. However, the empirical mean 
redshift of the UV dropout galaxies of Steidel et al. (1999) 
lies at z — 3.04 ±0.24 and here the above argument reverses, 
because there is then only a very limited redshift range in 
which a Lyman limit system can intervene in the U band. 
Thus at z ~ 3 there is the possibility that 80% of galaxies 
do not show UV dropout. Taking all the galaxies with simi- 
lar (r — i) cc d '■ {b — r) cc d colours to the UV dropout galaxies 
i.e. those with —0.1 < (b — r) cc d < 1.6 and (r — i) cc d < 0.2, 
we find ~ 220 new candidate z ~ 3 galaxies. Spectroscopy 
of these candidates are required to determine how many of 
them are actually at high redshift. 

The WHDF {u - b) ccd : (6 - r) ccd diagram in Fig. g(b) 
shows the same two populations as in (r — i) cc d '■ (b — r) cc d 
although slightly less distinctively. The plume of early types 
lies at (b — r) cc< j > 1.6 with the spirals now occupying a 
broad strip below them. Again the E/SO model seems to fit 
the E/SO plume toz~ 0.5 with the H/K break exiting the b 
band at z ~ 0.25 and then the evolved spectrum shortwards 
of 2500A entering the u band at z ^ 0.5 which takes the 
(u — b) cc d colour blueward until z ~ 0.8, when, with both 
u and b on the steep but constant slope of the UV 'upturn' 
spectrum, a constant (u — b) cc d = —0.4 colour is maintained 
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Figure 21. The K-corrections compared to the evolutionary tracks plotted in both (r — i) ccld : (b — r) cc( j and (u — b) ccd : (b — r) r ,.,i 
planes. Solid lines - E/SO K-correction; short dash - E/SO qo = 0.05 evolution; dots - Spiral K-correction; long dash - Spiral qo = 0.05 
evolution. Boxes indicate 0.5z steps for < z < 3.5. Small numerals indicate redshift. 




Figure 22. (a) the (r — i) cc d : (b — r) co( j colour-colour diagram for the WHDF data limited at b ccd < 25.75 compared with the predicted 
contours for the qo = 0.5 dE model. The contours are equally spaced in log surface density. (b) as (a) but for the qo = 0.05 model, (c) 
the data points shown without the contours for clarity - note the prominent 'hook' of spirals and the 'finger' of early-types. 



until the spectrum reddens again when Lyman a absorption 
enters the u passband at z 1.5. The strong bluewards 
shift in (6 — r) cc d between 0.7 ^ z ^ 1.5 is caused both 
by the H/K break being redshifted through r and the UV 
'upturn' spectrum at A < 2500A being redshifted through 
b. The spiral (u — b) ccd : (b — r) ccd tracks again first move 



bluer in [u — b) cc d and redder in (b — r) ccc i due to H/K 
moving through the b band at < z ^ 0.3. Then (6 — r) cc d 
becomes bluer at z <; 0.5 when the decreasing spectrum 
shortwards of the H/K break is redshifted through r until 
the UV 'upturn' reaches ratz = 1.5. Both u and b lie on the 
steeply increasing UV 'upturn' spectrum in the range 0.3 ^ 
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z Si 1.5 which maintains the (u — 6) cc< j roughly constant at 
(u — b) cc d ~ —0.8 before the spectrum reddens as dust and 
then Lyman a increasingly affects the u band at z <; 1.5. 

Also shown on Fig. pCj(b) (enclosed by the dotted lines) 
is the the area of the (u — b) cc d : (b — r) cc d colour-colour plane 
corresponding to the UV-dropout selection criteria of Steidel 
et al. (1999), as discussed in Section 6.3. The locations of our 
final sample of 43 u-band dropouts is indicated, with distinct 
symbols indicating the original 32 satisfying the criteria of 
Steidel et al., and those lost and gained from invoking the 
(r — i) cc d selection. Note that, as the dropouts are chosen 
from an r-limited sample, not all have corresponding data 
points on this {(-limited plot. It can be seen that there is 
good agreement with the location of z fa 3 galaxies implied 
by our models, for both early-types and spirals. 

As shown by Fig. [H], in (u — b) cc d : (b — r) C cd the effects 
of evolution are again bigger for the early-type model than 
the late-type model. At 0.5 z 1.8 the (6 — r) cc d colours 
for early-type evolution model become bluer than the K- 
corrections and at 0.5 ^ z 0.9 the (u—b) cc d colours for this 
model also become bluer than the K-correction. However, for 
2 > 0.9 in (u — b) cc d and z > 1.8 in (b — r) cc d the colours for 
the evolution model are redder than the K-correction and 
this trend is confirmed when intervening Lyman a absorp- 
tion enters both u and 6. For the spirals the K-corrections 
seem to give a similar result everywhere to the evolutionary 
model, evolution making the colours slightly bluer than the 
K-correction at low redshift and dust making the evolution- 
ary model slightly redder at high redshift. We note that the 
K-corrections of Pence (1976) are similar to the Bruzual & 
Chariot K-corrections for z S5 0.5 and, if anything, slightly 
bluer in (it — b) ccd at z i; 0.5. 

Thus in the redshift range 0.3 Si 2 Ss 1 the predicted 
(u — b) cc d colours for the spiral model are much bluer than is 
observed; the spiral model reaches (u — b) cc d = —1.0 whereas 
the data only reach (u — b) C cd = —0.5. The discussion of the 
Uccd and b cc d photometry zeropoints in section 3.2 indicates 
that errors in the colour zero-point of 0.1 mag or larger are 
unlikely. In any case we shall see below that the same effect 
is also seen in the independent HDF data. This discrepancy 
is surprising since the UV spectra of star-forming galaxies 
are only governed by two parameters, dust absorption and, 
more weakly, the IMF. We postpone to Section 7 further 
discussion of possible explanations of this result. 

6.5 HDF Colour-Colour Distributions 

An HDF-N colour-colour comparison with similar models 
has previously been briefly discussed by Metcalfe et al. 
(1996). We now add in the HDF-S data and look at the 
combined (F606 - F8U) vega : (F450 - F606) vega and 
(F300 - F450)„ e9<l : (F450 - F606) V ea a diagrams at both 
a bright and a faint limit in Figs. E3|. The most promi- 
nent feature remains the line of galaxies moving redward 
in (F300 - F4:50) vega while remaining blue in (F450 - 
F606)„e 9 a. These are the UV dropouts galaxies expected 
to be at z ii 2 as discussed above. In general, the tracks 
appear broadly consistent with the models and as such 
provide strong support for the approach used here. The 
(F606 - F8U) vega : (F450 - F606)„e B a tracks again seem 
good fits; the Sbc track in particular gives an excellent 
fit to many of the bluer galaxies, although when shifted 



to the rest colours of Scd and Sdm's gives slightly too 
blue a colour in (F450 — F606) vega . In the bright data the 
plume of early-type galaxies is not seen as clearly as in the 
WHDF in either (F300 - F450) vega : (F450 - F606)„ e9a 
or (F606 - F814)„ esa : (F450 - F606) vega but the models 
again give a reasonable representation of the redder galaxies. 
Again, as for the WHDF (u — b) cc d, the Sbc model predicts 
too blue (F300 — F450)„ e9tl colours for spirals at z — 1 and 
the arbitrary shifting of the Sbc curve to the rest colours of 
Scd and Sdm galaxies leaves their tracks not fitting well in 
this redshift range either. At fainter limits, it can be seen 
that galaxies predicted to lie at z ~ 2 dominate the data. 

Next, in Fig. ||(a),(b) we plot the HDF-N and HDF- 
S galaxies with spectroscopic redshifts on the (F300 — 
F450)„ e9 „ : (F450 - F606)„ esa and (F606 - F814)„ eg a : 
(F450 — F606)„ e9 a diagrams. The points are colour-coded 
according to the galaxy redshift range and the model tracks 
are similarly colour-coded. It can be seen there is good gen- 
eral agreement between the predicted and observed redshift 
ranges, again firmly supporting the idea that these simple 
models do an excellent job of describing the colours of faint 
galaxies. 

Fig. [25] shows how the HDF UV dropouts populate 
the correct part of the HDF (F606 - F814)„ e9a : (F450 - 
F606) vega diagram for z i; 2 galaxies. Again this indicates 
that ground-based (r — i) : (6 — r) colour-colour combina- 
tions may be effectively used to discriminate z <; 2 galaxies 
without the need of a UV band. 

Fig. ^ compares the predicted number density of galax- 
ies at each point in the (F606 - F814)„ e9a : (F450 - 
F606)i, e9 a colour-colour diagrams with the data at F606 < 
28. It can be seen that the qo=0.05 and qo=0.5 dE evolu- 
tionary models both give a good fit with the density peak- 
ing in both the models and the data at the colours of z ~ 2 
spirals. The fact that the spiral models fit so well through- 
out the range < z Ss 3 is again a strong argument that 
the evolutionary history of the SFR is well described by 
the r=9Gyr exponentially increasing SFR used here. The 
early-type galaxies are not so clearly seen in this HDF data 
compared to the equivalent WHDF data shown in Fig. 

6.6 Sizes 

It has been claimed that the angular sizes of faint HDF 
galaxies are small compared to local galaxies (Roche et al. 
1998 and references therein), suggesting there has been size 
as well as luminosity evolution. We have examined this ques- 
tion in the light of our simple evolutionary models by adapt- 
ing the simulation procedures used to replicate the WHT 
data in section 3.6 to produce mock HDF frames based on 
our models. We are then able to apply our data-analysis 
procedures to these frames and compare the angular sizes 
measured for the artificial galaxies with the real HDF data. 
As before (see Paper III), we assume the Freeman (1970) 
law to relate exponential disk size to absolute magnitude 
and the diameter-magnitude relations of Sandage & Perele- 
muter (1990) for r 1/4 bulges. Fig. || shows the results for 
our qo — 0.05 model. The solid line represents the mean 
isophotal radius (at /ZF450 = 28.5) of our artificial galaxies, 
whilst the dots show the same quantity for the F450-band 
HDF-N data. Although small differences exist, the agree- 
ment between model and data is remarkably good. It should 



© 2000 RAS, MNRAS 000, 



26 N. Metcalfe, et al. 




Figure 23. Minuted (F606-F814)„ escl : (F450 - F606)„ estl and (F300 - F450) vega ■ (F450-F606)„ ega colour-colour diagrams for the 
HDF-N and HDF-S data, split into bright and faint magnitude ranges. The evolutionary tracks predicted from our Bruzual & Chariot 
models are shown out to z = 4 for the five galaxy types. Filled squares show 2 = 0, open squares with numerals indicate 0.5z steps. 



be noted that the apparent isophotal angular radii for the 
faintest galaxies (~ 0.2") correspond to metric radii of only 
~ 2 kpc (Ho = 50 kms _1 Mpc _1 ). This is much smaller 
than local L* galaxies, but demonstrates both the effect of 
isophotal dimming and the fact that, at F450 vega ~ 29, in 
the models we are looking down the luminosity function to 
galaxies with low intrinsic luminosities and small sizes. 

We find similar results for the qo = 0.5 dE model, if 
we assume the "disappearing dwarf" population follows the 
same Sandage & Perelemuter law as local bulge dominated 
galaxies. 

6.7 Star Formation History 

Our evolutionary models imply that the SFR density rises to 
beyond z = 2, according to the r=9Gyr Bruzual & Chariot 



model. This is at odds with the conclusions of Madau et 
al. (1996) who used UV fluxes as an estimate of the SFR 
and, on the basis of the same HDF-N data, concluded that 
the SFR peaks at z ~ 1 and decreases rapidly thereafter. 
Following Shanks et al. (1998) we first discuss the reasons 
for the differences seen between our SFR history and the 
original SFR estimates of Madau et al. (1996). 

The approach of Madau et al. (1996), following Cowie et 
al. (1988) and Lilly et al.(1995,1996), uses the fact that the 
galaxy UV flux density is proportional to the SFR density to 
obtain estimates of the SFR history of the Universe. They 
use the CFRS spiral evolving LF of Lilly et al. (1995) to 
determine the 2800A flux density in the range < z Ss 1 
and the HDF F300W and F450W dropouts to obtain the 
1620A flux density in the ranges 2 Ss z Ss 3.5 and 3.5 ^ 
z 4.5. Their conclusion of a peak in the SFR at z w 1 is 
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(F606-F814) v „. (F300-F450) raes 

Figure 24. (a) (F606 - F814)„ e9a : (F450 - F606)„e 9 a and (b) (F300 - F450)„ ega : (F450 - F606)„ estl colour-colour diagrams for the 
HDF-N and HDF-S data showing the location of galaxies with known redshifts (coloured dots, with redshifts as indicated in the legend) 
compared with the model tracks out to z = 4. The tracks are colour-coded according to redshift range in the same fashion as for the 
galaxies. Arrows indicate limits to colours for galaxies not detected in one band, assuming a magnitude at the 3ct limit for the missing 
objects. 



Table 14. UV luminosity and SFR densities measured in the HDF and WHDF (H = 50fcms _1 Mpc' 1 ). LF 
indicates luminosity function incompleteness corrected. Ab = 0.3 indicates dust corrected. V e ff indicates 
that the effective volumes of Steidel et al. (1999) have been assumed, rather than the boxcar volume, V. In 
the WHDF case, factors are V/V e // for qo=0.5,0.05 respectively. 
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(go = 0.5, LF, 
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WHDF B< 26 (1 <, z < 2) 
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WHDF R< 25.5 (1 £ z £ 2) 


1.0 x 10 26 f 
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-1.55 


-1.87 


WHDF U dropouts(2.5 <i z & 3.5) 


3.0 x 10 25 


-2.55 


-2.80 


-1.94 


-2.18 


-1.61 


-1.85 


WHDF " (V e// , x2.38,3.15) 


7.1 x 10 25 


-2.17 


-2.30 


-1.56 


-1.64 


-1.23 


-1.31 


HDF-N U dropouts(2.0 <, z <; 3.5) 


2.1 x 10 26 


-1.70 


-1.95 


-1.40 


-1.73 


-1.07 


-1.40 


HDF-N " (Veff, xl.13) 


2.4 x 10 26 


-1.65 


-1.90 


-1.35 


-1.68 


-1.02 


-1.35 


HDF-S U dropouts(2.0 <, z < 3.5) 


2.5 x 10 26 


-1.62 
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-1.32 


-1.65 


-0.99 
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HDF-S " (Veff, xl.13) 


2.8 x 10 26 
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WHDF B dropout(3.5 <, z <J 4.5) 


2.2 x 10 25 


-2.68 


-2.94 


-1.98 


-1.92 


-1.65 


-1.59 


WHDF " (Veff, x 1.94,2.68) 


4.2 x 10 25 


-2.39 


-2.51 


-1.69 


-1.49 


-1.36 


-1.16 


HDF-N B dropout(3.5 <, z <, 4.5) 


8.1 x 10 25 


-2.11 


-2.36 


-1.89 


-2.17 


-1.56 


-1.84 


HDF-N " (V eff , xl.19) 


9.6 x 10 25 


-2.03 


-2.28 


-1.81 


-2.09 


-1.48 


-1.76 


HDF-S B dropout(3.5 <, z <, 4.5) 


7.6 x 10 25 


-2.14 


-2.33 


-1.92 


-2.14 


-1.59 


-1.81 


HDF-S " (V e f f , xl.19) 


9.0 x 10 25 


-2.06 


-2.25 


-1.84 


-2.06 


-1.51 


-1.73 



TP2800 



in contradiction with our t~9 Gyr exponentially increasing 
SFR for spirals, which we have shown fits a variety of faint 
galaxy count, colour and redshift data over a wide range of 
passbands. The SFR evolution rate of Madau et al. increased 
markedly faster than our exponential in the < z Si 1 range 
and then quickly decreased below our rate at z J> 1 (see fig. 
4 of Shanks et al. 1998 and also Fig. p3 below). 



At z = 0, Shanks et al. suggested that the reason our 
SFR density lies above that of Gallego et al. (1995) may 
be related to the bright B-band count normalisation issue, 
discussed above. Supported by the results of Glazebrook et 
al. (1995) and Driver et al. (1995), we explicitly ignore the 
low galaxy count at B ^ 17 mag as being possibly con- 
taminated by the effects of local large-scale structure. Fur- 
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Figure 25. The location of (F300 - F450)„ e9a drop-out galaxies in the HDF (F606 - F814)„ esa : (F450 - F606)„ e9 a plane. Circled 
dots indicate galaxies with no measured (F300 — F450) V ega- By comparing with Fig. 0(a) it can be seen that these drop-out galaxies 
have colours very close to those predicted by the models for galaxies with z > 2 in this colour plane. 



ther supporting evidence comes from the B < 19.5 mag 2dF 
galaxy redshift survey results of e.g. Colless (1999) where the 
luminosity function at low z appears to move in the density 
rather than the luminosity direction, as would be expected 
on the basis of the large-scale structure hypothesis. 

In considering possible reasons first for the differences in 
the 0.2 <, z <y 1 range, fig. 5 of Shanks et al. showed that our 
exponential models reasonably fit the rest B-band luminos- 
ity functions determined from the CFRS data. However, the 
worst fit of our evolving luminosity function to the CFRS 
data was in the 0.2 < z < 0.5 spiral bin and this proba- 
bly accounts for most of the discrepancy between our model 
SFR density and the CFRS data in this range. Further un- 
certainty is introduced because the 2800A flux density in 
this bin has to be extrapolated from the observed B-band 
data (Lilly et al. 1996) . Other measures of the SFR from the 
Oil measurements of Tresse & Maddox (1998) and Treyer et 
al. (1998) at z~0.2 agree much better with our exponential 
model at low redshift. 

In the 2^2^ 3.5 range, the reasons that the SFR 



density estimate of Madau et al. lies a factor of ~3 below 
our model in fig. 4 of Shanks et. al. (1998) are also clear. 
First, if we repeat the Madau et al. measurement in the 
HDF-N, then we obtain a factor of 1.56x bigger luminosity 
density at 1620A for z J; 2 UV drop-out galaxies than these 
authors, mostly due to the fact that we measure brighter 
magnitudes for individual galaxies (see Fig. ^). Second, if 
Madau et al. were to include our Ab(z=0) = 0.3 mag dust 
absorption for spirals, then this would effectively increase 
their 1620A luminosity density by another factor of 2.12. 
Combined, these two effects give a factor of 3.3 x overall 
which explains the above discrepancy of approximately a 
factor of 3. In the 3.5 z 4.5 bin the same two effects 
apply. However, the discrepancy here with our model is so 
large that we shall suggest below that they may in fact be 
in genuine disagreement. 

In Fig. ^(a,b) we illustrate some of these points by 
comparing various observational estimates of the SFR den- 
sity with our evolutionary models for both the qo=0.05 
(r=9Gyr, t ga ; aa; y=16Gyr, h 5 o = 1) and the q =0.5+dE 
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Figure 26. (a) the (F606-F814)„ e9a : (F450-F606)„ e ga colour-colour diagram for the HDF-N and HDF-S data, limited at F606„ e9a < 
28 compared with the predicted density contours of the go = 0.5 dE models. The contours are equally spaced in log surface density, (b) 
as (a) but for the qo = 0.05 model. 



(r=9Gyr, t ga i axy =12.7GyT, hso = 1) models. In (a) the dis- 
continuity arises from the contribution of the 'disappearing 
dwarfs with their constant SFR at z>l. The z — normali- 
sation used in both cases was obtained by first estimating the 
local luminosity density in the F300W passband from the 
Sbc/Scd/Sdm Schechter parameters in Table [l| This gave 
P3000A = 9.8 x 10 25 hw erg s" 1 Hz" 1 Mpc -3 . We then con- 
verted to SFR density using the factors given by Madau et 
al. (1996) (corrected downwards by a factor of 1.25 according 
to the extrapolation to 3000A from 1500-2800A and then 
downwards again by the factor 1.6 for the reasons suggested 
by Connolly et al. 1997) to obtain p, = 0.00735M©/i~ 3 yr" 1 
Mpc -3 . Finally we corrected this SFR(z) normalisation for 
our inclusion of Ab = 0.3 mag dust absorption. Assuming 
our 1/A dust model this amounts to an upwards correction 
by a factor of 1.54. 

All observational points based on P2800.4 are first cor- 
rected downwards by the factor of 1.6 recommended by 
Connolly et al (1997). They are then corrected using the 
Ab = 0.3 mag dust absorption as used in our model. Fur- 
ther, we have LF corrected each point either by taking the 
published correction or, where these are absent, as in the 
case of the HDF-N data of Madau et al. (1996), by run- 
ning the qo=0.05/ qo=0.5+dE evolutionary models in each 
redshift range and integrating the galaxy luminosity density 
to convergence. We also apply these LF corrections to our 
own data. Thus if our model and the data are fully consis- 
tent, then they should agree in this Figure. As well as the 
points of Gallego et al. and the CFRS points of Lilly et al. 
(1996) (both dust/LF corrected), we have also plotted the 



dust/LF corrected HDF-N points of Connolly et al. (1997) 
in the 1 Ss z 2 range. 

Also shown in Fig. ^ are our new WHDF estimates 
of the SFR density at 1 ^ z 2. These were obtained 
by integrating the luminosity density of galaxies with (b — 
r) C cd < 0.6 and (r — i) C cd > 0.4 which are approximately 
where the 1 ^ 2 2 galaxies lie in Fig. The 343 galaxies 
with b cc d < 26 and the 267 galaxies with r cc d < 25.5 give 
the values for the UV luminosity density and SFR densities 
shown in Table [w]; these are somewhat lower than given by 
Connolly et al. particularly in the qo=0.05 case. However, 
the colour cuts we have used are relatively crude and this 
may be responsible for any difference with the HDF results 
and our own data. 

In the redshift range 2.5 ^ z ^ 3.5, Fig. ^ shows 
the HDF-N UV dropout data of Madau et al. (1996), cor- 
rected for our A_g = 0.3 mag dust absorption and our 
LF incompleteness corrections. Also shown is the Keck 
UV dropout galaxy data of Steidel et al. (1999), with our 
dust (As = 0.3 mag with 1/A giving A1500 = 0.9 mag) 
and LF incompleteness corrections. With their own dust 
(E(B - V) C aizetu=0.15 giving A1500 = 1-7 mag) and LF 
incompleteness corrections, Steidel et al. (1999) find their 
z — 3 SFR is a factor of 3.4 higher than shown in Fig. ^ 
(a). Fig. ^ also gives new estimates of the SFR density at 
2.5 <; z <, 3.5 from WHDF u-band dropouts to r ccd <25.25 
together with the estimates from the data of Steidel et al. 
(1999). Here we have taken the 43 UV dropout candidates 
we identified in Section 6.4 to give the value of P1500 shown 
in Table [l4|. Again using the conversion factors stated in the 
paper of Madau et al. this gives the SFR densities shown 



© 2000 RAS, MNRAS 000, 



30 N. Metcalfe, et al. 



-0.5 



0.8 



0.6 



o 
.a 
a 

a 



0.4 



0.2 




' - -rv. i 



35 



26 



27 



28 



39 



F450 vsaa Magnitude 



Figure 27. The isophotal radii versus F450 ve ga magnitude re- 
lation for the HDF-N data (dots) compared with the prediction 
for the mean radii from a simulation of an HDF frame (solid line) 
based on our qo = 0.05 evolutionary model. 



in Table [w| The dust/LF corrected estimate is plotted in 
Fig. We have corrected our WHDF data upwards by 
a factor of 2.38 in the case of qo = 0.5, 3.15 in the case 
of qo = 0.05, to use the effective volumes (including their 
colour incompleteness/deblending correction) of Steidel et 
al. (1999). The observed WHDF SFR density estimate is 
close to the observed value from Steidel et al. (1999), when 
our corrections are consistently used. We similarly have plot- 
ted our estimates of the SFR at 2 Ss z Si 3.5 from the HDF- 
N and HDF-S UV dropouts using the same criteria as used 
by Madau et al. (1996). All HDF SFR's at z = 2.75 have 
been corrected upwards by a factor of 1.19 to account for 
colour incompleteness/blending as estimated by Steidel et 
al (1999). The HDF points agree well amongst themselves 
and with the model and are close to the WHDF and Steidel 
et al. points which have bigger LF incompleteness correc- 
tions. 

Further, we have plotted estimates of SFR density at 
3.5 Ss z Si 4.5 using F450W dropout galaxies from the 
HDF-N analysis by Madau et al. (1996), from the HDF- 
N and HDF-S analyses by ourselves using identical selection 
criteria, and from the WHDF, using the criteria of Steidel et 
al. (1999), transformed into our magnitude systems. These 
are all corrected using our dust and LF incompleteness cor- 
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Figure 28. (a) Star formation rate versus redshift in the case 
qo=0.5. The solid line represents our r = 9Gyr spiral evolution- 
ary model, normalised to include our dust absorption correction, 
with the contribution to the SFR density from the dE popula- 
tion added in at z> 1. This model is compared to data points 
from Lilly et al. (1996), Madau et al. (1996). Steidel et al. (1999), 
Gallego et al. (1995), Connolly et al. (1997) and from our esti- 
mates from the WHDF and the HDF. All data are corrected for 
LF incompleteness and Ab = 0.3 mag dust absorption. Where 
appropriate the SFR densities use the effective volumes of Steidel 
et al. (1999) (see Table [jj). (b) As for (a) for the case q =0.05. 
This time the solid line simply represents the t = 9Gyr spiral 
evolutionary model. 



rections (see Table |l4|). The WHDF data is also corrected 
upwards by the factor 1.94 in the case of qo = 0.5, 2.68 in 
the case of qo = 0.05, to use the effective volume of Steidel 
et al. (1999). All HDF data has also been corrected by a 
factor of 1.19 to use the effective volume recommended by 
Steidel et al. (1999). Also plotted are the SFR densities from 
Keck B dropout galaxies from Steidel et al. (1999); these are 
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based on 49 3.5 Si z ;S 4.5 galaxies with spectroscopic con- 
firmation. Again the points shown assume our dust and LF 
corrections. Using their own LF and dust absorptions would 
result in a factor of 1.7 increase in the SFR at z — 4 over the 
value shown in Fig. |28| (a). Although Steidel et al. include 

1.7 magnitudes of dust absorption at 1500A as opposed 
to our 0.9 magnitudes, their use of the effectively gray ex- 
tinction law of Calzetti (1997) means that they include 1.2 
magnitudes of dust at 2800A compared to our 0.5 magni- 
tudes, so the effect of their higher dust absorption is mainly 
on the normalisation of the SFR:z plot rather than its form 
(compare their fig. 9 to Fig. || (a)). In Fig. ||(a),(b) the 
WHDF points are close to both models at z = 4 but the Stei- 
del et al. point is higher than the qo=0.05 model at z = 4 
due to a large LF correction. However, in all cases the HDF 
points at z — 4 are lower, even though we have applied both 
the dust and LF incompleteness corrections consistently. We 
shall return to discuss this discrepancy between the HST and 
ground-based SFR estimates at z = 4 in Section 6.8. 

We conclude that, within the various errors and cor- 
rection factors, there is reasonable self-consistency with the 
fully corrected SFR density data and our model assuming an 
exponentially increasing SFR with look-back time, at least 
to z ~ 3. At 3.5 Ss 2 $j 4.5, the ground-based data of Steidel 
et al. and ourselves is generally in much better agreement 
with our exponential model than the HDF data. Steidel et 
al. (1999) have suggested that this may be because the HDF- 
N is too small to be representative. However, we note that 
the B dropout data from the HDF-S show good agreement 
with the HDF-N SFR density. Therefore the picture at high 
redshift based on the integrated measure of UV luminosity 
and SFR density remains unclear. We note that because of 
the large size of the LF incompleteness and dust corrections, 
even where our model and the SFR density data agree, this 
only really checks the internal consistency of our model and 
the SFR data. In Sect. 6.8 we shall make a more direct test 
of our models than afforded by Fig. 

Probes of the SFR at high z using SCUBA sub- 
millimetre observations (Smail et al. 1997, Hughes et al. 
1998, Barger et al. 1998) have also suggested that dust may 
re-radiate absorbed UV luminosity and thus provide an ex- 
planation for the sub-mm counts. Busswell & Shanks (2000) 
have shown that dust re-radiated UV flux in our evolving 
spirals, even though the B-band absorption is as low as 
0.3 magnitudes, can also make a significant contribution to 
the steep sub-mm counts. Thus our spiral evolution models 
which include dust may be enough in themselves to explain 
the sub-mm counts, at least at faint, ~ lmjy fluxes. 

6.8 Counts of UV and B dropout galaxies and the 
epoch of galaxy formation 

In order to clarify the status of our model predictions at 
2.5 z 4.5 we now consider Figures ^ and |30| where we 
have plotted the number-magnitude counts of UV dropout 
(as defined in section 6.4) and B dropout (see section 6.7) 
galaxies and compared them to our model predictions. The 
count predictions for the qo=0.05, qo=0.5(+dE) and A evo- 
lutionary models are shown and their shapes are similar to 
the galaxy luminosity functions at high redshift. The up- 
turn at M* in the qo=0.5 case is due to the extra contri- 
bution from the disappearing dwarf population at high red- 



Table 15. WHDF Galaxy counts for UV dropouts 2.5 < z < 
3.5 (from section 6.4). These are the raw data and have to be 
multiplied by a factor of 1.3 to account for colour incompleteness- 
blending after Steidel et al. (1999). 



Magnitude 

(r C cd) 

23.0-24.0 
24.0-24.5 
24.5-25.0 
25.0-25.25 



(per frame) (arcmin 2 ) 



22 
12 



0.01 ± 0.01 
0.16 ± 0.06 
0.45 ± 0.10 
0.49 ±0.14 



Table 16. HDF-N & HDF-S F606 vega galaxy counts for UV 
dropouts 2 < z < 3.5 (from Metcalfe et al. 1996). These are the 
raw data and have to be multiplied by a factor of 1.19 to account 
for colour incompleteness/blending after Steidel et al. (1999). 



Magnitude 




N 9 ai 


Magnitude 




N 9 ai 


(F606 vega ) 


(N) 


(arcmin 2 ) 


(F606„ e9a ) 


(N) 


(arcmin 2 ) 


HDF-S 






HDF-N 






23.0-23.5 


3 


0.56±0.19 


23.0-23.5 


1 


0.19±0.19 


23.5-24.0 





0.00 


23.5-24.0 





0.00 


24.0-24.5 


4 


0.75±0.38 


24.0-24.5 


4 


0.75±0.38 


24.5-25.0 


11 


2.63±0.70 


24.5-25.0 


8 


1.50±0.53 


25.0-25.5 


19 


3.56±0.82 


25.0-25.5 


21 


4.50±0.92 


25.5-26.0 


39 


7.31±1.17 


25.5-26.0 


35 


6.56±1.11 


26.0-26.5 


19 


3.56±0.82 


26.0-26.5 


28 


5.25±0.99 


26.5-27.0 


4 


0.75±0.38 


26.5-27.0 


7 


1.31±0.50 



shift. Taking Fig. |29|(a) first, which applies to the counts of 
ground-based UV dropouts in the redshift range 2.5 z 
3.5, plotted as a function of r cct j magnitude, we see that, 
although the new WHDF count (Table 09) is slightly lower 
than the Keck data from Steidel et al. (1999), overall the 
data agree well with the model predictions at bright mag- 
nitudes. Similarly, in Fig. ||(b), the HDF-N UV dropout 
data from Madau et al. (1996) and our own HDF-N and 
HDF-S analyses in the 2^2^ 3.5 range agree well with 
the models at fainter F606 vega magnitudes. In Fig. ^ we 
have adopted the average colour incompleteness/blending 
factors suggested by Steidel et al. (1999) which means that 
both the Steidel et al. and the WHDF number counts are 
increased by a factor of 1.3. We conclude that the luminosity 
function of galaxies in the redshift range 2 Sd z Si 3.5 ap- 
pears reasonably consistent with the local galaxy luminosity 
function of Sbc/Scd/Sdm galaxies (which dominate at high 
redshifts), as modified by the effects of r=9Gyr evolution 
plus dust. In the particular case of the qo=0.5 model, the 
data seem more consistent with a model which includes the 
dE population rather than the same model without the dE 
population (also shown in Fig. ^) . Thus there is strong ev- 
idence that the space density of galaxies at low redshift and 
z ~ 3 may be very similar, whatever the underlying model, 
and this may be a stronger conclusion than any drawn from 
Fig. ^ above which frequently have to rely on large LF in- 
completeness corrections at high redshift to compare SFR 
densities with low redshift measurements. Fig. |29| suggests 
that not only may the galaxy luminosity density be constant 
in the range 1 z 3 but the space density and luminosity 
function of galaxies may also be reasonably unchanging in 
this redshift range. It also suggests that the Lyman-break 
galaxies detected by Steidel et al. may be interpreted as the 
evolved bright tail of the local Sbc/Scd/Sdm galaxy lumi- 
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Figure 29. (a) The r cc[ j number count of 2.5 < z < 3.5 ground- 
based UV dropout galaxies, from Steidel et al. (1999) and from 
the WHDF, compared to our evolutionary model predictions with 
qo=0.05, 0.5 and f2A=0.7. In the qn=0.5 case the contribution of 
the dE population can be seen at r cc[ j > 25. A colour incom- 
pleteness/blending factor of 1.3 has been applied to the WHDF 
data, the same factor used to correct the data of Steidel et al. 
(1999). (b) The F606„ e9a number count of 2 < z < 3.5 HDF 
UV dropout galaxies compared to these same models. A colour 
incompleteness/blending factor of 1.19 has been applied to all the 
HDF data as suggested by Steidel et al. (1999). In all cases the 
agreement with the models is good. 



nosity function. These conclusions would not be altered by 
assuming the Calzetti (1997) dust law with E(B-V) = 0.15 
since the differential extinction between 1500 A at z = 3 and 
the B band at z = is similar to that found in our 1/A dust 
model with Ab — 0.3 mag. 

We now turn to Fig. ^ where the counts of B-band 
dropout galaxies from Steidel et al. (1999), Madau et al. 
(1996) in the HDF-N and ourselves in HDF-N and HDF-S 
and in the WHDF in the range 3.5 z ^ 4.5 are plotted 



Table 17. i-band differential galaxy counts for b- band dropouts 
on the WHDF. These are the raw data and have to be multiplied 
by a factor of 1.3 to account for colour incompleteness/blending 
after Steidel et al. (1999). 



Magnitude 

(iced) 

23.0-23.5 
23.5-24.0 
24.0-25.5 
24.5-25.0 



(per frame) (arcmin 2 ) 




3 
5 
10 



0.00 

0.061 ±0.04 
0.10 ±0.05 
0.20 ±0.06 



Table 18. Differential galaxy counts for FAbOW — band dropouts 
from the HDF-N and HDF-S fields. These are the raw data and 
have to be multiplied by a factor of 1.19 to account for colour 
incompleteness/blending after Steidel et al. (1999). 



Magnitude 






Magnitude 




N 9 ai 


(F814„ e9a ) 


(N) 


(arcmin 2 ) 


(F814„ es a) 


(N) 


(arcmin 2 


HDF-S 






HDF-N 






24.0-24.5 





0.00 


24.0-24.5 


3 
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as a function of i cc d magnitude. Again, the data is com- 
pared to the same four evolutionary models as in Fig. |29[ 
Although the three basic models still fit the brighter data 
of Steidel et al. (1999), they overpredict the fainter WHDF 
and HDF data by about a factor of 5. At z ~ 4, the data 
seem more similar to the fourth model shown, which is the 
qo=0.5 model without the extra dE population. Thus the 
overall conclusion seems to be that at z « 4 the observed 
galaxy luminosity function appears to drop significantly be- 
low our basic models at absolute magnitudes fainter than 
M*. 

One interpretation would be that the assumption of a 
constant space density of galaxies in our models had failed 
by z ~ 4 and that we were possibly detecting the epoch 
of galaxy formation. It might also be suggested in the con- 
text of the qn=0.5 model that the results might imply that 
the bright galaxies were in place even as early as z ~ 4 
and the dE population appeared between 3 & z 4 and 
then dimmed to become the dwarf galaxy population possi- 
bly detected at b < 25 in the WHDF (r — i) cc d : (b — r) ccd 
colour-colour diagram. This may be similar to the 'downsiz- 
ing' scenario claimed by other observers at lower redshifts 
(e.g. Cowie et al. 1996). 



7 DISCUSSION OF MODEL SUCCESSES AND 
FAILURES. 

Our basic conclusion is that the galaxy counts in the optical 
bands can be well fit by simple Bruzual & Chariot models 
with an exponentially increasing rate of star-formation to 
z ~ 3. In the low qo case, this type of model can fit the data 
to B ~ 27 whereas in the high qo case the model only fits to 
B ~ 25, just beyond the current limit for galaxy redshift sur- 
veys. In this case a model with either a disappearing dwarf 
dE population or an increased luminosity function slope at 
high redshifts must be invoked to increase the counts at faint 
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Figure 30. The i cc d number count of 3.5 < z < 4.5 B-band 
dropout galaxies from the ground-based data of Steidel et al. 
(1999) and from the WHDF and also from the HDF-N data of 
Madau et al. (1996) and the HDF-N and HDF-S data from this 
paper. A colour incompleteness/blending factor of 1.3 has been 
applied to the ground-based data and a factor of 1.19 to the 
HDF data. These data are compared to our evolutionary model 
predictions with qn=0.05, 0.5 and Q^=0.7. In the qn=0.5 case the 
contribution of the dE population can be seen at i cc d > 25 and 
the prediction for the qn=0.5 model without the dE population 
is also shown. The observed faint galaxy LF appears to fall below 
our model predictions by a factor of «5, indicating that while 
bright L* galaxies may have been in place at z fs4, the faint 
galaxy population may have formed between z 3 and z fs: 4. 



magnitudes. We have also shown that the models produce 
a magnitude-size relation, which for either value of qo, is 
consistent with the data. 

We have also presented new results from the WHDF 
and the HDF-N and HDF-S for the SFR density of galaxies 
over the redshift range 1 2 3.5 and the results from this 
approach are broadly consistent with both previous data and 
the above results if effects such as internal dust absorption 
in spirals are included in the models. However, other system- 
atic effects may still dominate this route to the SFR history. 
Also, with this approach the poor fit to the faint counts of 
the qo=0.5 model using the local galaxy luminosity function 
would be missed. 

At high redshifts, the number of galaxies predicted by 
the simple models at z ~ 3 agrees well with the number 
of UV dropout galaxies detected at R ~ 25 by ourselves 
and by Steidel et al. (1999). Thus the PLE models seem to 
continue to fit the data as far as z«3. However, at z « 4, the 
PLE models start to overpredict the faintest data and this 
may suggest that dwarf galaxies may have formed at z » 4. 
We regard the agreement between the spiral models and 
the bright, Lyman break galaxy data as another indication 
that the spiral model is broadly correct in the range 1 Si 
z SS 3 and therefore that the star-formation rate follows our 
exponential model throughout this redshift range. 

We have seen that the agreement between these models 



and the WHDF and HDF galaxy colours is particularly good 
in (R — I) : (B — R), where two faint galaxy populations are 
clearly seen corresponding to early- and late-type galaxies. 
The features seen in the data match those predicted by the 
models so well that this forms a strong argument that these 
models must be a reasonably accurate description of the 
evolving spectra of faint galaxies. Even in (U — B) : (B — R) 
they at least provide a good qualitative description of the 
data. 

The outstanding issue for the early-type galaxy models 
is the need to use a dwarf-dominated IMF (with a 0.5Mq 
cut-off) for early-types in order to reduce the number of z > 
1 galaxies predicted in K < 19 redshift surveys. The effect 
of using this IMF on the predicted optical colours is small. 
A slightly more extended initial period of star-formation is 
required with our dwarf dominated IMF(r=2.5Gyr) than 
with Salpeter (r=lGyr) to match the move towards bluer 
(B — R) colours at z ~ 0.5 seen in the (R — I) : (B — R) 
diagram. However, we note that the Salpeter/r=lGyr model 
would then also overpredict the numbers of z > 1 galaxies 
detected in optical (B < 24) redshift surveys too. Thus, over 
a wide range of wavelengths, the case against the Salpeter 
plus r=lGyr model for early-types from the redshift surveys 
is strong. The alternative solution to the dwarf-dominated 
IMF is a merging model as considered by other authors (e.g. 
Kauffmann & Chariot 1998). It remains to be seen whether 
the merging model will produce as good a fit to the (R — I) : 
(B — R) data as the simple PLE model discussed here. 

For the spiral model, the outstanding issue is that in 
the redshift range 0.5 Si z Ss 1 the predicted (u ~ b) cc d 
colours are up to 0.5 magnitudes too blue compared to 
the data. Since the u and b bands peak close together at 
3800A and 4400A this implies quite a sizeable discrep- 
ancy in the redshifted spectral slope between 1500-2500A , 
where the Bruzual and Chariot model predicts approxi- 
mately f\ oc A -2 ' 5 and the data would prefer f\ oc A ' 5 . 
(u — b) cc d is more affected than (b — r) cc d because the small 
width of the u band means there is less smoothing and 
because the wide wavelength separation of b and r means 
that they simultaneously lie in the range 1500-2500A for a 
shorter redshift span (at z ~ 1.7). The (it — b) cc d discrepancy 
is puzzling because the UV colours of star-forming galaxies 
are governed by only two parameters, dust and the IMF. We 
already have As = 0.3 mag dust in the model and because 
the reddening vector runs parallel to the (u — b) cc d : (b— r) cc d 
track at z = 1, it would take a significant increase in the 
amount of spiral dust to improve the fit and this would spoil 
the n(b) and n(z) fits even at b cc d ~ 25. 

We have experimented with including the Galactic 
2200A dust absorption feature but although this tends to 
help while it moves through u, it then moves through b mak- 
ing (u — b) cc d sharply redder and this tends to make the 
agreement with the data even worse at z <; 0.7. We also 
tried using an x=3 IMF as for the early-types with r=9Gyr 
to redden the (u — b) cc d colours but although this helped 
redden (it — b) cc d in the 0.5 Ss z SS 1 range, (b — r) cc d also 
became too red for any agreement with the data even at 
low redshift. We conclude that the original parameters of 
As = 0.3 mag, 1/A, dust absorption and a Salpeter IMF re- 
mains closer to the (u — b) cc d : (b — r) cc d data than these al- 
ternatives. The fact that the K-corrections seem to demand 
a steep UV slope for galaxies locally might be thought to 
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indicate an explanation for the {u — b) cc d problem where the 
IMF might evolve with redshift. However, the UV spectra 
of local spirals has much scatter and it has long been known 
that only the bluest local galaxies (e.g. NGC4449) give spec- 
tra as steep as predicted by the Bruzual & Chariot models 
with an IMF that is a single power-law. 

G. Bruzual has suggested (priv. comm) that an upper 
mass cut-off of 5-6 Mq in the spiral IMF rather than the 125 
Mq which we conventionally use will flatten the spiral s.e.d. 
shortwards of 2000A and help take the colours closer to the 
data. He has supplied models which show this effect on the 
(u—b) cc d '■ (b—r)ccd diagram and it clearly helps the (u—b) cc d 
colours at z = 1, while continuing to produce reasonable 
colours for a spiral at z — 0. A 5-6 Mq cut corresponds 
to an IMF with no stars earlier than about B2-B3 which 
might be plausible, particularly if more massive stars are 
assumed to be preferentially cocooned in dust. A potential 
problem with this or any other explanation which effectively 
reduces the slope of the UV s.e.d. for spirals is that the 
K -Revolution correction at a given SFR may produce less 
bright galaxies at high redshift, causing our evolutionary 
model to underpredict the B counts, for example. Increasing 
the SFR at early times to account for this might not be 
possible because the galaxy would evolve to be too red for 
a spiral at z — 0. 

Flat UV spectra have also been observed in the spectra 
of local starburst galaxies. For example, Gonzalez Delgado et 
al. (1998) have obtained UV spectra of four local starbursts 
with the Hopkins Ultraviolet Telescope and the Goddard 
High Resolution Spectrograph on HST. These authors show 
that the spectral energy distribution follows fx cx A* 3 with 
= -0.84 ± 0.26. This is much flatter from the f3 = -2.6 
to —2.2 they expect on similar grounds to our argument 
above. They, however, believe that the prediction for (3 is 
so stable that any deviation can be immediately interpreted 
as being caused by dust absorption. They obtain an average 
value of As = 0.73 ± 0.11 for their four starburst galaxies 
using an LMC extinction law. They also obtain the relation 
E(B - V) = 0.11A/3 where A/3 is the difference between 
the observed and predicted slope. In our case A/3 ~3 and so 
Ab ~ 1.3. This is so large that the dust would redden the 
galaxy colour predicted for an Sbc spiral by 0.6 magnitudes 
and a prime motivation for adopting the r = 9Gyr model, 
that it gave a blue enough colour for a present day spiral, 
would then be lost. Gonzalez Delgado et al. also detect SilV 
and CIV absorption lines with P Cygni profiles which is a 
signature of O-type and early B-type stars. They suggest 
that this argues against a bright IMF cut-off earlier than 
B2. 

Steidel et al. (1999) have fitted reddenings to their 
Lyman-break galaxies and find values very similar to those 
we use here, finding an average E(B — V) = 0.15 for an 
average Ab ~ 0.5 mag using the Calzetti (1997) absorp- 
tion law. However, E(B — V) = 0.15 would only produce 
a spectral slope change of A/3 = 1.2, and this is less than 
half what is needed to redden spiral (u — b) cc d colours to 
what is observed in the WHDF. Pettini et al. (2000) have 
also detected P Cygni profiles in the spectra of the lensed 
Lyman-break galaxy MS 1512-cB58 which again may indi- 
cate the presence of O-stars and argue against the existence 
of an upper mass cut-off. In this galaxy the spectral index 
implies E(B — V) = 0.10 again close to what we have as- 
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Figure 31. (a) The effect on our Sbc evolutionary model 
(qn=0.05) for (u — b) cc[ j : (b — r) cc[ j of a dust reddening law 
which evolves from a Galactic law exhibiting a 2200A absorp- 
tion feature at low redshifts (z < 0.6) to a 1/A extinction law at 
high redshifts (z > 0.6). The evolution helps to redden the model 
z = 1 to agree better with the data, (b) As for (a) showing the 
effect of the evolving dust model on the model predictions for 
(r - i)ccd ■ (b - r) ccd . 



sumed in our spiral models. However, this dust absorption 
implies a far higher dust-gas ratio than seen in comparably 
low-metallicity galaxies such as the LMC and SMC. 

Our view is that it may be unreasonable to invoke dust 
absorption as high as E(B — V) = 0.3 or Ab = 1.3 mag 
to explain the (u — b) cc d colours of the z — 1 spirals in the 
WHDF. On the other hand it is clear that P-Cygni profiles 
of SilV and CIV are seen in both Lyman-break galaxies and 
nearby starbursts. These features are seen in stars as late as 
B0 and a 5-6 M Q cut would correspond to a cut at B2-B3. 
We are therefore forced to consider another model where the 
evolution of the dust extinction law is now invoked to explain 
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the red (w— b) cc d colours of the WHDF data. Noting that the 
2200 A feature helped redden (u — b) C cd as it passes through 
the u band, we suggest that the high metallicity galaxies 
found at low redshift may therefore show a 2200A feature 
as seen in the Milky Way. Then, since the predicted sharp 
feature caused by the 2200A feature passing through the b 
band is not observed, we hypothesise that at higher redshifts 
than z « 0.6 the extinction law evolves to look more like 
the extinction law seen in the low metallicity SMC, which 
is close to our 1/A law. In this way, the strength of the 
2200A feature produces increased reddening as the u (and 
b) bands move towards 2200 A than in (b — r) cc d at higher 
redshift when the 2200A feature has disappeared. We have 
tested this evolutionary model using the following simple 
evolutionary prescription to mimic this effect in the u band. 

A x = A B x 4400/A, A > 2500A 

A\ = Ab x (4400/A + (9886/A - 3.9544) x (1 - 
exp{-t exp /t))), A < 25001 

where t exp =8Gyx and t is the look back time. This ensures a 
smooth transition from the 2200A feature dominated law at 
low z to the 1/A law at higher redshifts. This dust evolution 
then produces a much redder (it — b) cc d at z — 1 while not 
disturbing unduly the fit in (r — i) cc d : (b — r) cc d at either 
this or higher redshift (see Fig. [jlj ). 



8 CONCLUSIONS 

We have produced u cc d,b cc d,r cc d and i cc d galaxy number 
counts and (u — b) cc d, (b — r) cc d and (r — i) cc d colours from 
one 7' x 7' field, on which we have the deepest ground-based 
B-band image yet taken. We also have deep H and if-band 
photometry on this field (Paper IV, McCracken et al. in 
prep). 

We have compared these data with our analysis of the 
HDF-N and HDF-S, produced using the same software. For 
the UV and blue-bands, although the HDF's are deeper for 
unresolved objects, in terms of surface brightness signal-to- 
noise our WHT image is ~ 0.4 mag deeper, and also covers 
5 times the area. The counts from the datasets agree well, 
after the ground-based data has been corrected for confu- 
sion. All are still rising at the faintest magnitudes, with a 
slope dlog(N)/dm ~ 0.25. The HDF F606W counts are over 
2 magnitudes deeper than the previous equivalent ground- 
based limits. They also show a still-rising count at the faint 
end with a slope dlog(N)/dm ~ 0.25. The ground-based U 
counts show an almost constant slope of dlog(N)/dm ~ 0.4 
to U ~ 25.5. However the HDF F300W counts flatten to a 
slope dlog(N)/dm ~ 0.15. 

The deepest counts are the HDF-S STIS data, which 
reach to STISab ~ 30, which is roughly equivalent to R ~ 
30. Our evolutionary models continue to fit to the limit of 
these deepest optical number count data. 

In terms of the comparison with evolutionary models, 
we summarise our conclusions as follows: 

1. We have fitted simple exponential star- formation his- 
tory models to our data, and found that low qo cosmologies 
give a good fit from U to K over a range of ~ 10 magnitudes. 
To enable high qo models to fit, an extra component must 
be added to the galaxy population at high (z <; 1) redshift. 
One way of doing this is by adding a 'disappearing dwarf 



dE population which is bright and blue at high redshift and 
faint and red at low redshift. 

2. At very bright magnitudes (B ^ 17) our count models 
fail and overpredict the counts. Since the steep count at 
B <, 17 seems to be reproduced at the same level at R, I, 
H and K, which is unlikely for an evolutionary explanation, 
this represents strong evidence that the high normalisation 
we have adopted for the counts is substantially correct. 

3. The models provide a good fit to the distributions of 
colours in (6 — r) cct j and (F450 — F606) vega and in (r — i) cct j 
and (F606 — F81<i) vega , although in both (u — b) cc d and 
(F300 — -F450)„e 9 a the spiral models are too blue by 0.5 
magnitudes at z — 1. 

4. The same models produce elliptical and spiral tracks 
in the colour-colour planes with shapes very similar to the 
loci of the data. The fact that the fit of our models is so good, 
particularly in (r — i) cc d '■ (b—r) C cd for the ground-based data 
and (F 606 - F814)„ ega : (F 450 - F606) vega for the HDF, 
data lends strong support to the underlying assumption that 
the star-formation rate increases exponentially to beyond 
z = 2 (Metcalfe et al. 1996). It is clearly possible to use 
(R — I) : (B — R) to probe the star-formation rate directly 
over a wide redshift range 0.5 z & 3, emphasising the 
power of this colour-colour diagram as a diagnostic tool for 
galaxy evolution. 

5. In the WHDF (r — i) C cd '■ (b — r) cc d plane, we detect 
an excess of low redshift(z « 0.1), low luminosity early- 
type galaxies which may have parameters consistent with 
the 'disappearing dwarfs' invoked in our qo=0.5 evolutionary 
model. 

6. We have shown how the discrepancy with the z = 1 
SFR peak of Madau et al. (1996) comes from the differ- 
ent treatment of dust and photometry in the HDF. When 
we repeat the UV dropout SFR density measurement made 
in the 2 Ss z is 3.5 range by Madau et al. we find good 
self-consistency with the prediction of our exponentially in- 
creasing SFR (r=9Gyr) spiral model. 

7. We have made new measurements of the SFR den- 
sity in the WHDF and in the HDF-N and HDF-S. In the 
redshift range 1 <, z <, 2 the WHDF SFR density lies lower 
than the measurement of Connolly et al. (1997) in the HDF. 
These measurements bracket the prediction of our r=9Gyr 
spiral model. In the 2 ^ z £ 3.5 range our observed SFR 
density based on UV dropouts, corrected for dust and LF 
incompleteness, and the ground-based estimates of Steidel 
et al. (1999) lie close to the HDF estimates and our model. 
At 3.5 z 4.5 there is a wider spread, with the brighter 
ground-based SFR estimates from the WHDF and Steidel 
et al. generally being higher and closer to the model than 
the HDF SFR estimates which are significantly lower than 
the model. 

8. We have detected a population of galaxies that oc- 
cupy the same part of the (r — i) C cd ■ (b — r) cc d plane as 
z ~ 3 UV dropout galaxies, but which do have measurable 
u-band fluxes. We have suggested that these may be candi- 
dates for high redshift galaxies that do not show UV dropout 
and that statistics of QSO absorption lines may not rule out 
such a population at z ~ 3. The increase in SFR density at 
z ~ 3 will depend on the fraction which prove to be at high 
redshift. 

9. We have created simulated CCD frames to check 
the apparent sizes of galaxies predicted by our evolutionary 
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models and found they are quite consistent with the sizes of 
galaxies measured on the HDF-N. 

10. For the early-type models, the main remaining prob- 
lem is that although they produce reasonable fits to the 
data at optical wavelengths, they predict too many high 
redshift (z 1) galaxies at K ~ f9 unless they have a 
dwarf-dominated stellar IMF. This reduces the evolutionary 
brightening at infra-red wavelengths whilst still allowing suf- 
ficient optical evolution. Dynamical merging may constitute 
an alternative explanation and this possibility has been ex- 
plored by other authors (eg Kauffmann & Chariot 1998). 

11. The main outstanding question for the spiral mod- 
els is why do they predict (u — b) cc d and (F300 — -F450)„ esa 
colours which are too blue in the range 0.5 ^ z & 1? We 
have suggested that this problem is unlikely to be solved 
by the addition of more dust to the models because the 
amount of dust needed is so great (Ab ~ 1.3 mag) that it 
would lead to the model predicting (B — R) colours that 
are ~0.5 magnitudes too red for az = spiral. One possi- 
bility is that the spiral IMF is truncated above 5-6 MQ. It 
is the subject of continuing research to see if such a model 
can produce enough evolution to continue to fit the steep 
B-band number-count slope. Another possible explanation, 
suggested by the dust laws in nearby galaxies with different 
metallicities, is that the dust extinction law evolves, with the 
2200A dust absorption feature appearing in high metallic- 
ity galaxies at z Si 0.5 and disappearing in lower metallicity 
galaxies at higher redshift. 

12. The evolutionary qo=0.05 and qo=0.5 models cor- 
rectly predict the sky density of Lyman break z w 3 galax- 
ies detected by Steidel et al. (1999) from Keck spectroscopy. 
The models suggest that these galaxies may be spiral galax- 
ies which have undergone an exponentially increasing SFR 
with look-back time. At z ~ 4 these models generally over- 
predict the space density of sub-L* galaxies at faint limits 
in the HDF-N and HDF-S and WHDF by a factor of 5, sug- 
gesting that we may have reached an epoch of significant 
galaxy formation. The z ~ 4 luminosity function matches 
best our qo=0.5 model with the dE population removed, 
suggesting that one interpretation might be that bright L* 
galaxies were already in place at z = 4 but many fainter 
galaxies formed at 3 ?S z £ 4. These then quickly dimmed 
to become the red, dwarf population possibly detected in the 
WHDF (r - i) ccd : (b - r) ccd data at M s = -15 at z » 0.1. 
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